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ABSTRACT 


Two  irregular  wave  tests  were  conducted  over  a  permeable  fine  sand  beach 
in  equilibrium,  recording  time  series  of  the  free  surface  and  horizontal  velocity  at  17 
cross-shore  locations.  The  first  part  of  this  study  supplements  the  results  of  Herrman 
et  al.  (1997),  obtained  using  a  1:16  impermeable  smooth  slope.  The  exponential 
gamma  distribution,  with  measured  mean,  standard  deviation,  and  skewness,  is 
shown  to  be  capable  of  describing  the  probability  distributions  of  the  above  data  in 
a  unified  manner  over  the  beach  profiles,  in  mean  water  depth  h  >  \  cm.  In  the 
swash  zone,  the  skewness,  s,  approaches  2  and  the  root-mean-square  wave  height 
Hims  approaches  a/S  h.  Cross-shore  variations  of  the  mean  and  standard  deviation 
of  mid-depth  horizontal  velocity  are  adequately  predicted  in  the  shoaling  and  surf 
zones  using  linear  long-wave  theory  locally.  In  addition  to  cross-shore  data,  two 
velocity  profiles  were  measured  and  horizontal  velocity  statistics  were  shown  to  be 
fairly  constant  over  depth.  An  attempt  is  also  made  to  synchronize  free  surface  and 
velocity  time  series  to  obtain  information  about  incident  and  reflected  wave  spectra 
in  shallow  water.  Total  wave  spectra  are  examined  to  check  for  the  existence  of  long 
period  standing  waves  over  the  profiles.  Finally,  sheet  flow  and  energetics-based 
sediment  transport  models  are  used  with  the  measured  horizontal  velocities  to  try 
to  predict  the  equilibrium  profiles  obtained  in  this  experiment.  The  energetics  model 
is  shown  to  predict  these  profiles  more  accurately  than  the  sheet  flow  model.  In  the 
milder  wave  conditions  of  the  second  test,  the  energetics  model  predicts  equilibrium 


within  experimental  uncertainty  of  1  cm/hr,  but  both  models  significantly 
overpredict  bar  movement  during  the  more  intense  wave  breaking  of  the  first  test. 
Their  shortcomings  are  examined  so  as  to  suggest  possible  improvements  in  the 


models. 
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Chapter  1 


INTRODUCTION 


1.1  The  Swash  Zone 

The  swash  zone  on  a  beach  is  located  between  the  surf  and  backshore  zones. 
It  is  the  area  of  the  beach  that  is  alternately  submerged  and  exposed  by  wave  action. 
A  predictive  model  for  swash  dynamics  is  required  in  order  to  model  the  erosion  and 
recovery  of  natural  and  artificial  beaches.  Swash  data  has  been  taken  on  natural 
beaches,  thereby  providing  a  qualitative  understanding  of  swash  dynamics  [e.g., 
Guza  and  Thornton  (1982);  Holman  and  Sallenger  (1985);  Douglas  (1992)].  Yet  it  is 
unlikely  that  simple  empirical  relationships  will  be  developed  from  this  data  because 
of  the  limitless  possible  combinations  of  different  beach  profiles  and  incident  wave 
characteristics.  The  time-dependent  numerical  model,  based  on  the  finite-amplitude 
shallow- water  equations  including  bottom  friction  (Kobayashi  and  Wurjanto  1992), 
has  been  shown  to  be  capable  of  predicting  the  swash  characteristics  on  a  natural 
beach  within  errors  of  about  20%  (Raubenheimer  et  al.  1995).  However,  this  shallow- 
water  model  cannot  be  extended  to  intermediate  water  depth.  Furthermore,  the 
time-dependent  numerical  model  requires  significant  computation  time  to  resolve 
the  breaking  wave  profiles  that  vary  rapidly  in  time  and  space.  As  an  alternative  to 
the  above  procedures,  it  may  be  worthwhile  to  examine  possible  statistical  methods 
for  predicting  swash  dynamics. 


1.2  Previous  Statistical  Work 

Time-averaged  models  in  which  random  waves  are  represented  by  the  root- 
mean-square  wave  height  [Battjes  and  Janssen  (1978);  Thornton  and  Guza  (1983)] 
or  expressed  as  the  superposition  of  regular  waves  [Dally  and  Dean  (1986);  Mase 
and  Kobayashi  (1991)]  are  more  efficient  computationally  than  time-dependent  mod¬ 
els.  However,  they  sacrifice  detailed  temporal  information  such  as  the  skewness  of 
the  wave  profile,  which  was  found  to  be  important  in  predicting  cross-shore  sedi¬ 
ment  transport  [e.g.^  Guza  and  Thornton  (1985)].  Furthermore,  the  existing  time- 
averaged  models  may  considerably  underestimate  the  wave  setup  and  root-mean- 
square  wave  height  in  the  swash  zone  (Cox  et  al.  1994). 

In  estimating  nearshore  sediment  transport  rates,  Guza  and  Thornton  (1985) 
computed  the  various  moments  of  the  fluid  velocity  field  measured  on  a  gently  slop¬ 
ing  beach.  They  assumed  that  time  averaging  and  probabilistic  averaging  of  a  ran¬ 
dom  variable  were  equivalent.  Under  this  assumption,  comparisons  of  the  computed 
moments  with  monochromatic  and  linear  (Gaussian)  random  wave  models  indicated 
major  shortcomings  in  the  monochromatic  representation  of  random  waves.  The 
Gaussian  model,  with  given  mean  and  standard  deviation,  was  shown  to  predict  the 
even  moments  fairly  accurately  but  to  be  incapable  of  predicting  the  odd  moments 
associated  with  nonlinearities.  For  relatively  deep  water,  attempts  were  made  to 
predict  the  skewed  distribution  of  the  free  surface  elevation  using  higher-order  ex¬ 
pansions  [e.g.,  Longuet- Higgins  (1963);  Huang  and  Long  (1980);  Tayfun  (1980)].  In 
shallow  water,  Bitner  (1980)  compared  a  Gram-Chalier  expansion  with  field  data. 
Still,  this  higher-order  distribution,  with  measured  skewness  and  kurtosis,  predicted 
the  measured  free  surface  probability  density  functions  only  slightly  better  than  the 
Gaussian  distribution. 
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1.3  Previous  Work  With  Beach  Profiles  and  Sediment  Transport 

At  the  present  time,  there  are  no  models  capable  of  accurately  predicting 
the  development  and  final  shape  of  terraced  and  barred  equilibrium  beach  profiles. 
Profiles  in  the  form  of  d  =  with  d  =  still  water  depth,  y  =  seaward  dis¬ 

tance  from  the  shoreline,  and  A  =  sediment  scale  parameter,  have  been  shown  to 
represent  typical  profiles  of  natural  beaches  [e.g.^  Dean  (1991)].  Perturbations  of 
these  equilibrium  profiles  in  the  form  of  barred  and  terraced  beaches  have  also  been 
observed.  However,  the  processes  involved  in  the  creation  of  these  perturbed  beach 
profiles  are  not  presently  well  understood.  Trowbridge  and  Young  (1989)  developed 
a  sheet  flow  sand  transport  model  based  on  the  bottom  shear  stress  and  showed 
that  this  model  could  explain  the  measured  onshore  movement  of  a  nearshore  bar 
in  Duck,  North  Carolina  during  the  mild  wave  conditions  between  February  and 
August,  1982.  Thornton  et  al.  (1996)  used  an  energetics-based  sediment  transport 
model  to  explain  the  offshore  movement  of  a  bar  6n  the  same  beach  during  storms 
in  October,  1990.  Yet  no  existing  model  can  predict  both  onshore  and  offshore  bar 
migrations  satisfactorily. 

1.4  Objectives  and  Overview 

In  recent  years  the  probabilistic  approach  has  been  neglected  in  favor  of  time- 
dependent  numerical  modeling,  probably  for  lack  of  a  reasonably  simple  skewed 
distribution.  The  following  study  is  in  part  an  attempt  to  test  the  accuracy  of  one 
such  distribution  and  to  provide  experimental  data  for  the  construction  of  a  time- 
averaged  model.  In  this  respect,  it  is  a  continuation  of  work  begun  by  Herrman  et 
al.  (1997). 

In  their  work  and  in  this  research,  the  exponential  gamma  distribution  [e.g., 
Gran  (1992)]  with  given  mean,  standard  deviation,  and  skewness  is  adopted  to  de¬ 
scribe  the  probability  density  function  of  the  free  surface  elevation  from  outside  the 
surf  zone  through  the  swash  zone.  This  distribution  allows  skewness,  s,  to  range 
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between  s  =  0  and  5  =  2.  The  skewness  determines  the  shape  of  the  distribution, 
which  becomes  Gaussian  for  s  =  0  and  exponential  for  s  =  2.  The  Gaussian  dis¬ 
tribution  may  be  appropriate  in  relatively  deep  water,  far  seaward  of  the  breaker 
zone,  where  local  nonlinearity  may  be  negligible.  The  exponential  distribution  may 
be  a  good  approximation  in  the  swash  zone,  where  the  lower  limit  of  the  free  sur¬ 
face  elevation  is  determined  by  the  elevation  of  the  bottom.  As  there  has  been 
no  theoretical  justification  for  the  adopted  exponential  gamma  distribution,  Her- 
rman  conducted  experiments  on  a  uniform  1:16  slope  to  evaluate  its  applicability,  as 
well  as  to  examine  the  cross-shore  variations  of  the  mean,  standard  deviation,  and 
skewness  of  the  free  surface  elevation  from  outside  the  surf  zone  to  the  swash  zone. 
The  fixed  1:16  slope  eliminated  the  uncertainty  of  the  bottom  elevation  but  did  not 
represent  typical  concave  beach  profiles. 

This  study  conducts  similar  tests  with  a  concave  fine-sand  beach,  on  one 
barred  and  one  terraced  equilibrium  profile.  It  compares  the  predictions  of  the 
exponential  gamma  distribution  to  surface  elevation  and  cross-shore  velocity  data 
and  examines  the  cross-shore  variation  of  free  surface  and  velocity  statistics  across 
these  irregular  profiles.  Velocity  statistics  over  the  depth  are  also  examined  at  two 
cross-shore  locations  to  determine  the  degree  of  vertical  variation  of  the  velocity  in 
the  region  above  the  bottom  boundary  layer  and  below  the  wave  trough  level. 

The  second  major  objective  of  this  study  is  to  explore  the  shortcomings  of 
existing  sediment  transport  theories  in  predicting  perturbed  equilibrium  profiles. 
The  sheet  flow  and  energetics-based  models  discussed  in  section  1.3  are  applied  to 
the  equilibrium  profiles  obtained  in  the  two  tests  in  order  to  examine  the  reasons 
why  existing  sand  transport  models  cannot  predict  these  equilibrium  profiles  of 
no  net  cross-shore  sediment  transport.  The  measured  free  surface  elevations  and 
velocities  are  analyzed  to  obtain  the  free  surface  and  velocity  statistics  included  in 
these  models.  The  sheet  flow  model  is  only  capable  of  predicting  onshore  sediment 


transport  and  does  not  include  the  effects  of  undertow.  The  energetics-based  model 
has  been  shown  to  predict  only  offshore  bar  movement  and  growth  due  to  strong 
undertow  and  longshore  currents.  Suggestions  will  be  made  for  the  development  of 
an  improved  model  to  explain  the  observed  equilibrium  profiles. 

The  following  chapter  introduces  the  exponential  gamma  distribution  and 
corresponding  predictions  for  free  surface  parameters  in  the  swash  zone.  Chapter  3 
provides  an  overview  of  the  laboratory  experiments,  including  initial  testing  of  the 
fine  beach  sand  and  establishment  and  measurement  of  the  equilibrium  profiles.  It 
also  discusses  the  wave  gauge  and  velocimeter  probe  locations  and  provides  a  basic 
description  of  the  two  wave  tests.  Chapters  4  and  5  present  an  analysis  of  measured 
free  surface  elevations  and  velocities  respectively,  obtaining  the  probability  distri¬ 
butions  of  these  variables  and  comparing  velocity  results  to  the  predictions  of  linear 
wave  theory.  This  data  is  then  supplied  to  the  sheet  flow  and  energetics  sediment 
transport  models  in  Chapter  6.  The  deviations  of  these  model  predictions  from 
the  observed  zero  net  transport  may  provide  a  guide  to  developing  more  accurate 
predictions  of  cross-shore  sediment  transport. 


5 


Chapter  2 


EXPONENTIAL  GAMMA  DISTRIBUTION 

This  chapter  develops  the  theory  of  the  exponential  gamma  distribution  and 
describes  some  of  its  implications  for  basic  wave  parameters  in  the  lower  swash  zone, 
which  remains  wet  during  wave  uprush  and  downrush.  It  is  assumed  in  this  study 
that  the  time  series  of  the  free  surface  elevation,  rj,  above  the  still  water  level  (SWL) 
is  a  stochastic  process  satisfying  the  conditions  of  stationality  and  ergodicity  [e.^., 
Goda  (1985)].  The  time  and  probabilistic  averaging  of  the  stochastic  variable  can 
then  be  regarded  to  be  equivalent,  and  they  will  thus  both  be  denoted  by  an  overbar 
in  the  following.  The  probability  density  function  of  ?/  is  assumed  to  be  expressed 
by  the  exponential  gamma  distribution  \e.g.,  Gran  (1992)].  The  mean,  standard 
deviation,  and  skewness  coefficient  of  7/  are  denoted  by  g,  a,  and  s,  respectively. 
It  is  assumed  that  0  <  a  <  2  on  the  basis  of  available  field  and  laboratory  data 
[Bitner  (1980);  Huang  and  Long  (1980);  Goda  (1985);  Mase  and  Kobayashi  (1991); 
Raubenheimer  et  al.  (1995)]. 

2.1  Form  of  the  Exponential  Gamma  Distribution 

The  normalized  free  surface  elevation,  ?;*,  is  defined  as  g*  =  {f}  —  v)Ict  and 
thus  has  mean  and  standard  deviation  equal  to  zero  and  unity,  respectively.  The 
skewness,  a,  of  77*  is  the  same  as  the  skewness  of  rj.  The  exponential  gamma  distri¬ 
bution  for  77*  can  be  shown  to  be  expressed  as 


with 


y  =  rj^  -  (2.2) 

in  which  a  =  shape  parameter;  F  =  gamma  function;  V’  =  digamma  function;  and 
=  trigamma  function.  The  relationship  between  s  and  a  is  given  by 

s  =  -r{a)  [liV)]''-'  (2.3) 

in  which  =  tetragamma  function.  The  kurtosis,  jpf ,  of  77*  and  rj  is  given  by 

K  =  r'(«)  +  3  (2.3) 


in  which  ■0'"  =  pentagamma  function. 

The  gamma  and  related  functions  are  explained  in  Abramowitz  and  Stegun 
^  (1972)  and  tabulated  by  Gran  (1992).  The  values  of  r(a),  0(a),  0'(a),  0"(a),  and 

0'"(a)  for  given  a  can  be  found  using  Mathematica  (Wolfram  1991).  The  value  of  a 
for  given  s  can  be  obtained  by  solving  (2.3)  using  an  iteration  method  (Gran  1992). 
•  Figure  2.1  shows  the  variations  of  a  and  (TF  — 3)  with  respect  to  s.  For  s  <  0.4,  a  can 

be  described  approximately  by  a  =  (0.5  +  s“^)  (Gran  1992).  The  shape  parameter, 
a,  decreases  from  infinity  at  s  =  0  to  zero  at  s  =  2.  The  kurtosis,  K ,  increases  from 
_  Fl  =  3  at  s  =  0  to  /F  =  9  at  s  =  2. 


2.2  Limits  of  the  Exponential  Gamma  Distribution 

In  deep  water,  skewness  approaches  zero  and  equation  (2.1)  reduces  to  the 
Gaussian  distribution: 

In  the  wet  lower  swash  zone,  skewness  approaches  two  and  equation  (2.1)  becomes 
the  exponential  distribution; 

f(r)^,)  =  for  s  =  2  (2-6) 
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This  exponential  distribution  is  limited  to  the  range  ?/*  =  [{rj  —  ^)/cr]  >  —1.  This 
lower  limit  may  be  interpreted  in  terms  of  the  instantaneous  water  depth,  h  = 
[t]  —  Zb)^  where  zi  =  a,  bottom  elevation  which  is  taken  to  be  positive  above  SWL. 
Here,  the  bottom  elevation  is  assumed  to  be  in  equilibrium  and  to  remain  unchanged 
in  each  test.  It  is  required  that  h  >  0  and  hence  that  f]  >  zi,.  The  lower  limit  of  r] 
for  s  =  2  may  be  assumed  to  be  imposed  by  resulting  in  {fj  —  cr)  =  Zb-  Since  the 
mean  water  depth  h  is  given  by  ^  =  (^  —  Zb),  this  intuitive  assumption  requires  that 
the  mean  depth  approach  the  standard  deviation  in  the  swash  zone: 

h  =  (T  for  s  =  2  (2-7) 


2.3  Trends  of  Available  Data 

As  in  Herrman  et  al.  (1997),  one  objective  of  this  study  is  to  investigate 
whether  equations  (2.6)  and  (2.7)  hold  in  the  swash  zone.  Previous  research  has 
been  inconclusive  and  sometimes  contradictory.  Earlier  surf  zone  laboratory  data 
[e.g.,  Mase  and  Kobayashi  (1991)]  indicated  a  landward  decrease  in  skewness,  s, 
toward  the  still  water  shoreline.  Raubenheimer  et  al.  (1995)  measured  the  free 
surface  oscillations  using  runup  wires  placed  parallel  to  a  natural  beach  face  and 
obtained  small,  negative  values  of  s  for  the  sea  swell  frequency  range.  However, 
the  data  of  Herrman  et  al.  (1997)  showed  that  the  skewness  did  indeed  increase 
landward  of  the  still  water  shoreline  and  approach  a  value  close  to  two  in  the  swash 
as  predicted  above.  This  study  will  provide  further  evidence  in  support  of  (2.6)  and 
(2.7).  Note  that,  as  in  Herrman  et  al.  (1997),  the  swash  zone  free  surface  elevation, 
77,  used  below  corresponds  to  that  measured  by  a  vertical  gauge  with  its  lower  limit 
on  the  fixed  beach  face. 

The  root-mean-square  wave  height,  Hrms,  is  defined  in  this  study  using  the 
standard  deviation,  a,  [Thornton  and  Guza  (1983);  Battjes  and  Stive  (1985)] 


If  equation  (2.7)  holds  in  the  swash  zone,  then  Hrmsih  =  y/S,  which  is  much  larger 
than  the  result  Hrmsih  —  0.42,  measured  in  the  inner  surf  zone  of  a  natural  beach  by 
Thornton  and  Guza  (1982).  On  the  other  hand,  the  free  surface  elevation  time  series 
measured  by  vertical  capacitance  wave  gauges  buried  partially  in  the  sand  beach  in 
the  SUPERTANK  Project  (Kriebel  1994)  indicated  that  Hmolv  ranges  between  3 
and  4  at  the  still  water  shoreline,  where  Hmo  =  zero-moment  wave  height,  given 
by  Hmo  =  4:(T.  Since  Zb  =  0  and  h  =  rj  at  the  still  water  shoreline,  this  observed 
range  corresponds  to  H„aslh  =  2.1  to  2.8.  This  suggests  that  equation  (2.7)  may  be 
reasonable  in  the  swash  zone  and  that  the  value  of  Hrms! h  may  increase  appreciably 
from  the  inner  surf  zone  to  the  swash  zone.  The  data  of  Herrman  et  al.  (1997)  on 
a  1:16  impermeable  slope  also  indicated  that  Hrrasjh  approached  a  value  near  3  in 
the  swash  zone.  This  study  will  provide  additional  support  for  these  results. 


Chapter  3 


EXPERIMENTAL  PROCEDURES 


Because  velocity  statistics  were  measured  in  multiple  positions  in  this  study 
using  a  single  velocimeter,  it  was  necessary  to  run  the  same  random  wave  set  repeat¬ 
edly  over  an  unchanging  equilibrium  profile,  repositioning  the  velocimeter  between 
runs.  For  this  reason,  it  was  important  both  to  achieve  an  equilibrium  beach  profile 
with  the  fine  sand  used  in  the  experiment  and  to  obtain  repeatable  sets  of  waves.  In 
order  to  predict  the  beach  profile  shape,  it  was  necessary  to  obtain  further  informa¬ 
tion  about  characteristics  of  the  sediment.  Special  techniques  were  required  for  the 
calibration  and  processing  of  swash  wave  gauges  to  ensure  accurate  measurement  of 
relatively  small  water  depth  above  the  sand  bottom.  This  chapter  will  review  these 
experimental  procedures,  summarize  their  results,  and  describe  the  experimental 
setup. 

3.1  Sediment  Characteristics 

Approximately  6  m^  of  fine  white  quarry  sand  was  used  in  the  construction  of 
an  initially  1:12  beach  slope.  To  anticipate  better  the  shape  of  the  equilibrium  pro¬ 
files  that  would  be  established  from  this  initial  slope  and  to  facilitate  later  analysis 
of  sediment  transport,  several  tests  were  conducted  to  determine  basic  characteris¬ 
tics  of  the  fine  sand.  These  tests  measured  density,  specific  gravity,  porosity,  grain 
size  distribution,  and  fall  velocity. 


Density,  specific  gravity,  and  porosity  were  estimated  using  standard  testing 
procedures  (A1  Khafaji  and  Andersland  1992).  Bulk  density  was  measured  by  de¬ 
termining  the  weight  of  a  known  volume  of  lightly  compacted  sand.  To  find  the 
specific  gravity  this  measured  quantity  of  sand  was  then  added  to  a  known  volume 
and  weight  of  water  in  a  flask,  while  maintaining  the  same  total  volume  in  the  flask. 
It  was  then  heated  for  ten  minutes  to  remove  trapped  air,  and  then  weighed.  The 
specific  gravity  of  the  sand  was  defined  as  the  ratio  of  the  immersed  sand  weight 
to  that  of  the  pure  water.  Sand  water  content  was  determined  to  be  insignificant. 
Using  specific  gravity  and  water  content  data,  porosity  was  calculated  as  the  vol¬ 
ume  of  voids  divided  by  the  gross  volume  of  the  sand.  The  bulk  density  of  the  sand 
was  found  to  be  p  =  1.7  kg/m^.  Specific  gravity  was  determined  to  be  s  =  2.66, 
indicating  somewhat  silty  sand.  Porosity  was  measured  as  rip  =  0.4. 

3.1.1  Grain  Size 

Grain  size  was  measured  using  a  sieve  test  procedure  with  seven  different 
sieve  sizes.  The  numerical  results  of  three  sieve  tests  are  are  shown  graphically  in 
Figure  3.1  and  summarized  in  Table  3.1.  Sample  sizes  with  total  mass  (My)  ranging 
from  206.1  g  to  499.2  g  were  agitated  for  15-25  minutes.  The  relative  steepness  of 
the  three  curves  in  Figure  3.1  indicates  that  the  sand  is  of  fairly  uniform  size  and 
thus  will  establish  an  equilibrium  profile  with  little  or  no  variation  of  sand  size  along 
the  wave  tank.  As  shown  in  the  figure,  mean  diameter,  dso,  was  determined  by  linear 
interpolation  of  the  portions  of  the  three  curves  between  sieve  sizes  0.150  mm  and 
0.212  mm.  In  each  case,  dso  was  found  to  be  0.18  mm. 


Percent  passing 


Table  3.1:  Summary  of  Sieve  Test  Data 


Sieve 

Size  (mm) 

Sant 

Mass 

M 

i  capture 

Percent 

of  Mt 

■4(1) 

Percent 

Passing 

Sand 

Mass 

(s) 

capture 

Percent 

of  Mx 

d(ii) 

Percent 

Passing 

Sand 

Mass 

(9) 

capturec 

Percent 

of  Mx 

1(111) 

Percent 

Passing 

0.250 

33.4 

6.87 

93.13 

12.9 

6.26 

93.74 

20.1 

4.03 

95.97 

8.9 

4.32 

89.42 

16.2 

3.25 

92.72 

MilIrfiB 

85.8 

41.63 

47.79 

221.3 

44.33 

48.39 

0.125 

35.1 

7.22 

36.90 

15.8 

3.17 

45.22 

117.9 

24.25 

■HM 

BililfiM 

39.9 

8.21 

■291 

10.92 

5.29 

(Pan) 

21.6 

4.44 

0.00 

10.9 

5.29 

0.00 

99.2 

11.22 

0.00 

TOTALS 

486. 2g 

100.00 

206.  Ig 

100.00 

- 

499.2g 

100.00 

- 

3.1.2  Fall  Velocity 

By  equating  the  drag  and  the  submerged  weight  of  a  spherical  particle  falling 
in  quiet  water,  one  obtains  the  relation 

IM  "  '1^ 

in  which  w  =  particle  terminal  fall  velocity,  g  =  gravitational  constant  =  9.81  m/s^, 
s  =  particle  specific  gravity,  d  =  particle  diameter,  and  Cd  =  drag  coefficient. 
Average  fall  velocity  for  the  fine  sand  was  calculated  by  substituting  s  =  2.66  and 
dso  =  0.18  mm  into  equation  (3.1),  along  with  an  empirical  value  of  Cd  for  a  sphere 
[e.g.,  Yalin  (1976)].  It  was  estimated  that  w  =  1.89  cm/s  for  the  spherical  particle. 

Fall  velocity  was  also  determined  experimentally.  A  few  grains  from  each  size 
group  in  the  first  sieve  test  were  allowed  to  fall  through  a  clear  plastic  cylinder  filled 
with  water,  and  their  motion  was  visually  timed  for  a  distance  of  one  meter.  Ten 
fall  times  were  recorded  for  each  size  group  and  then  averaged  to  obtain  a  mean  fall 
velocity  for  each  group.  Results  are  summarized  in  Table  3.2.  When  these  results 
were  averaged  with  weighting  proportional  to  each  group’s  percent  of  the  total  mass. 


Table  3.2:  Fall  Velocity  Measurement  Results 


Sieve 

Number 

Effective 

Diameter 

Percent 
of  Mfotal 

Measured 

w 

60 

0.27mm 

6.26% 

5.2  cm/s 

70 

0.230 

4.32 

3.0 

100 

0.178 

41.63 

2.1 

120 

0.137 

8.44 

1.6 

140 

0.115 

23.14 

1.3 

200 

0.089 

10.92^ 

1.0 

Pan 

— 

5.29 

0.67 

it  was  found  that  Wavg  =  1.9  cm/s,  practically  the  same  as  the  1.89  cm/s  predicted 
for  the  sphere  of  diameter  d^o  —  0.18  mm. 

3.2  Experimental  Setup 

The  experiments  were  conducted  in  a  wave  tank  that  was  30  m  long,  2.44  m 
wide,  and  1.5  m  high  with  a  constant  water  depth  of  61.0  cm.  Repeatable  irregular 
waves,  based  on  the  TMA  spectrum  (Bouws  et  al.  1985)  using  linear  wave  theory  and 
random  phases,  were  generated  with  a  piston-type  wave  paddle.  The  experimental 
setup  is  shown  in  Figure  3.2.  A  rock  beach  with  a  1:8  slope  was  located  at  the  far 
end  of  the  tank  to  absorb  waves  propagating  outside  of  the  experimental  area.  The 
previously  constructed  divider  wall  along  the  center  line  in  the  tank  was  again  used 
to  enclose  the  beach,  and  a  rear  wall  was  added  to  support  the  upper  beach  berm 
area.  As  mentioned  above,  the  fine  sand  beach  was  built  with  an  initial  slope  of 
1:12.  Ten  capacitance  wave  gauges  were  used  to  measure  the  temporal  variations  of 
free  surface  elevations.  An  Acoustic-Doppler  Velocimeter  (ADV)  (Kraus  et  al.  1994) 
was  employed  to  measure  the  temporal  variations  of  fluid  velocities.  The  sampling 
rate  for  the  free  surface  and  fluid  velocity  measurements  was  20  Hz.  Two  tests 
were  conducted  on  two  different  sand  profiles  after  equilibrium  had  been  established 
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Figure  3.2:  Wave  Tank  Experimental  Setup 
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under  the  specified  random  wave  action.  The  wave  paddle  motion,  controlled  by  a 
computer,  was  approximately  identical  for  all  the  runs  in  each  test. 


3.3  Establishing  Equilibrium  Profiles 

To  estimate  whether  acceptable  equilibrium  profiles  would  be  attained  for 
the  given  sand  beach  and  incident  waves,  use  was  made  of  the  simple  profile  shape 
equation  introduced  by  Bruun  (1954); 

d  =  (3.2) 


where  d  =  still  water  depth  to  profile  (m),  y  —  distance  offshore  from  still  water 
shoreline  (m),  and  A  =  sediment  scale  parameter  based  on  the  mean  diameter  of 
the  sand.  It  is  noted  that  d  is  used  in  this  study  for  both  sediment  diameter  and 
still  water  depth,  in  accordance  with  their  customary  notations.  For  sand  with  dso 
=  0.18mm  it  was  estimated  that  A  =  0.094  m^/^.(Dean  et  al.  1994). 

Initial  attempts  to  predict  onshore/offshore  sand  transport  in  this  study  made 
use  of  the  profile  parameter  for  random  waves  in  the  field  (Kraus  and  Mason  1993): 

p  =  (3.3) 

where  g  =  9.81  m/s^,  Hs,o  =  significant  wave  height  in  deep  water  (m),  w  =  sediment 
fall  velocity  =  0.019  m/s  in  this  experiment,  and  Tp  =  peak  period  of  incident 
random  waves.  For  values  of  P*  >  26,500,  data  generally  shows  offshore  transport 
and  the  creation  of  bars  in  the  field.  For  the  random  wave  sets  chosen  to  be  used 
in  tests  1  and  2,  the  incident  wave  spectral  peak  periods  were  Tp  =  2.8  s  and  1.6  s 
and  the  significant  wave  heights  were  Hs,o  =  0.21  m  and  0.18  m,  respectively.  For 
these  wave  conditions,  the  corresponding  profile  parameters  were  determined  to  be 
approximately  22,000  and  29,000.  This  implied  empirically  that  there  would  be  a 
tendency  for  onshore  transport  with  the  formation  of  a  berm  in  test  1  and  offshore 
transport  with  the  creation  of  a  bar  in  test  2. 


Before  each  of  the  two  tests  was  begun,  the  random  wave  set  preselected  for 
each  test  was  run  repeatedly  for  several  days,  and  the  beach  was  allowed  to  shape 
itself  into  a  stable  equilibrium  profile.  Testing  commenced  only  after  the  profile 
shape  was  observed  to  remain  nearly  constant  for  approximately  one  day  of  wave 
action.  A  comparison  of  the  two  final  equilibrium  profiles  for  tests  1  and  2  is  shown 
in  Figure  3.3,  along  with  a  plot  of  the  predicted  profile  shape.  While  the 


Figure  3.3:  Final  Equilibrium  Profiles  (-  -  Test  1;  Test  2)  and  d  =  Ay^^^ 

Ay2/3  profile  follows  the  general  trend  of  the  two  equilibrium  profiles,  it  does  not 
account  for  the  significant  perturbations  in  the  region  of  the  bars.  The  steeper  slope 
further  offshore  is  partly  a  result  of  the  initial  1:12  slope  which  is  steeper  than  typical 
offshore  slopes.  In  partial  agreement  with  the  predictions  of  the  profile  parameter, 
the  test  1  profile  stabilized  as  a  terrace,  with  only  minor  evidence  of  a  bar  and 
trough.  A  larger  and  more  visible  bar  and  trough  developed  on  the  test  2  profile, 
after  a  period  of  offshore  transport. 

3.4  Measurement  of  Equilibrium  Profiles 

During  each  of  the  two  tests  detailed  measurements  were  made  of  the  beach 
profile  using  a  manual  vernier  pointer  in  the  swash  zone  and  a  Panametrics  22DLHP 
ultrasonic  depth  gauge  in  deeper  water.  To  ensure  that  longshore  profile  variation 


was  negligible,  the  profile  was  measured  along  three  cross-shore  transects  separated 
by  28  cm  in  the  longshore  direction,  with  the  center  transect  running  down  the 
center  of  the  experimental  beach  area.  Each  transect  included  77  measurement 
points  spaced  at  16  cm.  For  test  1,  three-transect  profiles  were  taken  before  the 
test  began,  in  the  middle  of  the  test,  and  after  the  test  was  complete.  For  test  2, 
because  of  complications  with  the  ultrasonic  depth  gauge,  it  was  only  possible  to 
take  a  two-transect  profile  before  the  test  and  a  three-transect  profile  after  the  test. 

As  shown  in  Figures  3.4  -  3.8,  longshore  variation  in  the  raw  profile  data 
was  generally  small,  especially  in  the  swash  zone.  The  cross-shore  distance  in  these 
figures  corresponds  to  the  shoreward  distance  from  the  coordinate  origin  at  the 
location  of  wave  gauge  1,  as  shown  in  Figure  3.2.  The  greatest  variation  occurred 
slightly  offshore  of  the  swash  zone  and  amounted  to  about  10%  of  the  water  depth. 
In  both  tests  fairly  uniform  ripples  were  established  in  the  regions  offshore  of  the 
bar  and  along  the  terrace  between  the  bar  and  the  swash  zone.  These  ripples  were 
approximately  1-3  cm  in  height  and  10-15  cm  in  length  and  were  observed  to  shift 
their  position  slightly  during  runs.  The  representative  profile  for  each  measurement 
set  was  obtained  by  first  smoothing  the  three  transects  with  a  three-point  smoothing 
routine  and  then  averaging  the  resulting  profiles.  The  final  equilibrium  profile  for 
test  1  was  obtained  by  averaging  the  pre-test,  mid-test,  and  post-test  smoothed 
profiles.  The  final  equilibrium  profile  for  test  2  was  an  average  of  the  pre-test  and 
post-test  smoothed  profiles  only.  The  smoothed  profiles  for  each  date  are  shown 
in  Figures  3.9  and  3.10.  The  small  profile  shifts  that  occurred  over  several  hours 
of  testing  are  represented  by  the  differences  between  these  lines.  Such  insignificant 
change  after  several  hours  of  wave  attack  confirms  that  the  profiles  were  essentially 
in  equilibrium.  The  final  profiles  for  tests  1  and  2  were  plotted  earlier  in  Figure  3.3. 


Figure  3.4:  Unsmoothed  Profiles  at  Three  Transects,  3/7/97  (Pre-Test  1) 


Figure  3.5:  Unsmoothed  Profiles  at  Three  Transects,  3/11/97  (Mid- Test  1) 


Figure  3.6:  Unsmoothed  Profiles  at  Three  Transects,  3/13/97  (Post-Test  1) 


Figure  3.7:  Unsmoothed  Profiles  at  Two  Transects,  5/2/97  (Pre-Test  2) 
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e  3.8:  Unsmoothed  Profiles  at  Three  Transects,  5/30/97  (Post-Test  2) 


Figure  3.9:  Smoothed  Average  Profiles,  3/7/97,  3/11/97,  and  3/13/97  (Test  1) 


Figure  3.10:  Smoothed  Average  Profiles,  5/2/97,  5/30/97  (Test  2) 


3.5  Gauge  Measurement  Positions 

Each  test  consisted  of  21  runs  performed  for  the  free  surface  measurements 
along  the  center  line  of  the  beach  profile.  Velocity  measurements  were  made  at  a 
location  approximately  10  cm  to  one  side  of  the  center  line  (see  Figure  3.2).  The 
duration  of  each  run  was  400  s,  from  which  the  initial  transient  duration  of  75  s  was 
removed.  The  remaining  325  s  of  each  time  series  was  used  for  the  subsequent  data 
analyses.  The  runs  from  each  test  were  numbered  1  through  21. 

3.5.1  Wave  Gauge  Positions 

The  cross-shore  locations  of  all  wave  gauge  measurement  positions  are  shown 
in  Figure  3.11,  over  the  two  final  equilibrium  profiles.  Wave  gauges  1,  2,  and  3  were 
fixed  for  each  test  and  employed  to  separate  the  incident  and  reflected  waves  using 
the  procedure  described  by  Kobayashi  and  Raichle  (1994).  The  gauges  were  located 
at  the  water  depth  d  =  60.0  cm  below  SWL  for  both  tests.  In  accordance  with 
the  procedure  of  Herrman  tt  al.  (1997),  the  horizontal  coordinate  x  in  this  paper  is 
taken  to  be  positive  landward  with  a;  =  0  at  the  location  of  gauge  1.  For  both  tests, 
wave  gauges  2  and  3  were  located  at  a:  =  0.38  m  and  x  =  0.50  m,  respectively. 
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Figure  3.11:  Cross-Shore  Positions  of  Wave  Gauge  Measurements  Over  Final  Pro¬ 
files  (-  -  Test  1,  Test  2) 


As  indicated  in  Figure  3.2,  wave  gauges  4  and  5  were  attached  to  a  movable 
carriage  along  with  the  ADV.  They  were  repositioned  17  times  during  each  test. 
Measurements  made  by  these  wave  gauges  were  spaced  at  0.25  m  and  covered  a 
region  from  x  =  2.60  m  to  10.85  m,  with  corresponding  depth  range  d  =  15-58  cm 
for  test  1  and  d  =  17-57  cm  test  2.  The  ADV  and  gauge  5  were  at  the  same  cross¬ 
shore  position  during  each  run  in  order  to  allow  for  synchronized  measurements 
of  surface  elevation  and  cross-shore  velocity.  Gauges  6  through  10  were  partially 
buried  in  the  sand  and  were  moved  only  once,  between  runs  12  and  13  for  test  1  and 
between  runs  10  and  11  for  test  2.  The  two  positions  of  gauge  6  were  spaced  0.70  m 
apart,  while  the  swash  measurements  of  gauges  7-10  were  spaced  at  0.14  m.  Water 
depths  at  gauges  4-10  for  the  two  tests  are  listed  in  Tables  3.3  and  3.4.  Note  that 
the  water  depth  d  is  used  for  gauges  4-6,  whereas  the  bottom  elevation  zi,  =  —d  is 
given  for  gauges  7-10  to  avoid  implying  a  negative  water  depth  in  the  swash  zone, 
where  d  <  0.  The  bottom  elevations  at  gauges  7-10  were  measured  by  a  manual 
vernier  pointer  for  each  run  because  the  relatively  small  water  depth  at  these  gauge 
locations  was  significantly  affected  by  the  minor  bottom  elevation  change  during 


runs. 


Table  3.3:  Depth  of  Water  at  each  Gauge  Position  for  Free  Surface  Test  1 


Water  Depth,  d 

(cm) 

Bottom  Elevation,  zi,(cm) 

Run 

Gauge  4 

Gauge  5 

Gauge  6 

Gauge  7 

Gauge  8 

Gauge  9 

Gauge  10 

■n 

58 

58 

14 

-6.6 

-3.4 

2.1 

7.4 

■a 

57 

56 

14 

-6.7 

-3.4 

2.0 

7.4 

1-3 

56 

55 

14 

-6.6 

-3.4 

1.9 

7.3 

1-4 

54 

53 

14 

-6.5 

-3.4 

1.8 

7.2 

1-5 

51 

50 

14 

-6.6 

-3.4 

1.7 

7.1 

1-6 

49 

47 

14 

-6.7 

1.6 

7.0 

1-7 

45 

43 

14 

-6.8 

1.7 

7.1 

1-8 

40 

37 

14 

-6.9 

1.7 

7.1 

1-9 

34 

32 

14 

-6.9 

-3.4 

1.6 

1-10 

28 

24 

14 

-6.8 

-3.5 

1.5 

OH 

28 

24 

14 

-6.7 

-3.5 

1.4 

6.9 

ng 

24 

-6.8 

-3.6 

1.4 

6.8 

1-13 

HEHI 

21 

11 

-6.5 

-2.3 

4.0 

10.0 

mi 

22 

23 

11 

-6.5 

-2.4 

3.8 

9.9 

—El 

23 

23 

11 

-6.3 

-2.3 

3.6  1 

9.8 

■nfiH 

22 

21 

11 

-6.1 

-2.3 

3.6 

9.8 

mam 

21 

21 

11 

-6.2 

-2.3 

3.5 

9.6 

1-18 

21 

21 

11 

-6.1 

-2.4 

3.5 

9.5 

1-19 

21 

21 

11 

-5.9 

-2.5 

3.4 

9.4 

1^ 

11 

-6.3 

-2.4 

3.2 

9.3 

1-21 

17 

15 

11 

-6.6 

-2.6 

3.1 

9.2 

Table  3.4:  Depth  of  Water  at  each  Gauge  Position  for  Free  Surface  Test  2 


Water  Depth,  d 

(cm) 

Bottom  Elevation,  Zh(cm) 

Run 

Gauge  4 

Gauge  5 

Gauge  6 

Gauge  7 

Gauge  8 

Gauge  9 

Gauge  10 

HI 

57 

57 

16 

-8.8 

-4.9 

0.4 

4.9 

55 

54 

16 

-9.8 

-4.6 

0.5 

5.0 

53 

51 

16 

-7.7 

-4.9 

0.0 

4.5 

2-4 

50 

48 

16 

-7.3 

-5.2 

0.0 

4.6 

2-5 

46 

44 

16 

-9.7 

-4.7 

0.4 

5.0 

41 

16 

-9.5 

-4.8 

0.4 

5.1 

35 

^HIH 

16 

-9.5 

-4.9 

0.4 

28 

26 

16 

-8.0 

-4.7 

0.2 

4.7 

2-9 

24 

22 

16 

-7.9 

-4.9 

4.7 

2-10 

22 

24 

16 

-7.7 

-4.9 

4.5 

mi 

22 

24 

14 

-3.1 

2.0 

6.4 

lina 

22 

24 

14 

imEH 

2.0 

6.3 

2-13 

25 

25 

14 

-5.3 

-3.0 

1.1 

6.2 

25 

25 

14 

-5.6 

-3.0 

2.0 

6.2 

24 

23 

14 

-5.6 

-3.1 

2.0 

6.2 

23 

22 

14 

-5.5 

-3.1 

1.9 

6.2 

2-17 

21 

20 

14 

-3.1 

6.2 

2-18 

21 

20 

14 

-5.5 

-3.0 

1.9 

6.1 

mi 

21 

14 

-3.0 

1.9 

6.1 

mi 

19 

18 

14 

-5.7 

-3.1 

1.9 

6.1 

2-21 

18 

17 

14 

-5.8 

-3.1 

2.0 

6.1 

Wave  gauges  6  through  10  were  partially  buried  in  the  sand  beach,  as  il¬ 
lustrated  in  Figure  3.2.  Gauges  7  through  10  were  located  within  the  swash  zone, 
while  gauge  6  was  in  water  that  was  too  shallow  to  allow  operation  of  an  ordinary 
capacitance  gauge.  Figure  3.12  depicts  the  special  design  of  these  swash  gauges, 
intended  to  minimize  scour  around  the  point  of  contact  between  gauge  and  sand. 
The  gauge  ground  was  connected  to  a  flat  steel  plate  and  attached  to  the  divider 
wall  shown  in  Figure  3.2,  while  the  active  end  of  the  gauge  consisted  only  of  a  thin 
insulated  rod. 


Figure  3.12:  Swash  Wave  Gauge  Design 


3.5.2  Wave  Gauge  Calibration 

Wave  gauges  were  calibrated  before  and  in  the  middle  of  each  of  the  two 
tests.  Calibrations  were  conducted  by  raising  the  level  of  water  in  the  tank  until  the 
sand  at  gauge  10  was  at  least  5  cm  under  water,  then  gradually  draining  the  tank 
and  recording  gauge  voltage  readings  for  every  1  cm  change  in  the  water  surface. 
For  every  calibration  25  different  readings  were  recorded  for  each  of  the  ten  gauges. 
For  the  pre-test  calibrations  the  water  level  was  initially  raised  to  -f  13  cm  relative 
to  the  normal  still  water  level,  while  for  the  mid-test  calibrations  the  water  level 
was  initially  raised  to  -1-15  cm.  This  level  was  then  allowed  to  drop  by  24  cm  to 
-11  cm  (pre-test)  or  -9  cm  (mid-test).  The  higher  initial  water  level  for  the  mid¬ 
test  calibration  was  necessary  to  obtain  a  sufficient  number  of  gauge  10  calibration 
points  after  that  gauge  was  moved  to  its  highest  position. 


Tank  Water  Level  Relative  to  SWL  (cm) 


Figure  3.13:  Calibration  Points  and  Linear  Fit,  Gauge  1,  Pre-Test  1 


For  gauges  1-5  calibration  data  generally  followed  a  straight  line,  as  illus¬ 
trated  for  gauge  1  in  Figure  3.13.  For  these  gauges  the  conversion  factor  was  ob¬ 
tained  from  the  best  fit  slope  to  all  of  the  data  points.  This  slope  is  shown  in 
Figure  3.13  as  a  solid  line.  The  vertical  dashed  line  in  the  figure  represents  the 
elevation  of  the  normal  still  water  level  during  the  wave  tests.  For  the  buried  gauges 
6-10  the  sand  bottom  was  eventually  exposed  by  the  falling  water  level  in  the  tank 
during  the  calibration.  At  this  point  readings  from  each  gauge  became  more  erratic 
and  strayed  from  the  straight  line  slope,  as  shown  for  a  gauge  9  calibration  in  Fig¬ 
ure  3.14.  In  this  plot  the  voltage  readings  follow  the  linear  fit  quite  well  while  the 
water  level  is  above  the  sand.  However,  they  immediately  depart  from  the  straight 
line  once  the  falling  water  exposes  the  sand  (indicated  here  by  the  vertical  dashed 
line).  This  departure  is  partly  due  to  a  delayed  drop  of  the  water  table  in  the  sand 
as  compared  to  the  water  level  in  the  tank.  It  is  also  likely  that  gauge  readings 
are  distorted  in  the  partially  saturated  sand  because  of  the  different  capacitance  of 
the  material.  Because  of  such  irregular  effects  in  partially  saturated  sand,  it  was 
decided  to  obtain  conversion  factors  for  the  buried  gauges  by  using  only  the  slope 
of  the  linear  fit  to  the  data  above  the  sand  level.  Such  a  slope  is  shown  by  the  solid 
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Figure  3.14:  Calibration  Points  and  Linear  Fit,  Gauge  9,  Pre-Test  1 

line  in  Figure  3.14.  In  order  to  obtain  at  least  five  linear  points  above  the  sand  at 
the  highest  position  of  gauge  10,  Zb  =  10  cm,  it  was  necessary  to  raise  the  water  in 
the  tank  to  its  highest  level,  15  cm  above  SWL. 

Following  test  1,  it  was  determined  that  readings  for  gauge  7  were  erratic  in 
both  calibrations.  Significant  variation  had  occurred  even  before  the  falling  water 
level  had  exposed  the  sand  at  that  gauge.  Figure  3.15  shows  an  example  of  this 
irregular  behavior  between  the  two  vertical  dashed  lines,  which  here  represent  still 
water  level  and  sand  level,  respectively.  Because  such  irregularity  would  have  re¬ 
sulted  in  additional  errors  of  up  to  3  cm  in  free  surface  measurements,  it  was  decided 
to  discard  all  gauge  7  records  from  the  data  analyzed  in  test  1.  The  first  twelve  runs 
for  gauge  8  in  test  1  were  also  discarded  for  the  same  reason.  For  test  2,  no  such 
calibration  irregularities  were  observed.  However,  it  was  necessary  to  discard  one 
record  each  for  gauge  5  at  ®  =  10.35  m  and  gauge  6  at  a;  =  12.25  m  in  run  20  of  the 
second  test  because  they  exhibited  large  and  unrealistic  deviations  from  all  other 
wave  records  in  the  test. 
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Tank  Water  Level  Relative  to  SWL  (cm) 


Figure  3.15:  Calibration  Points  and  Linear  Fit,  Gauge  7,  Mid-Test  1 
3.5.3  Velocity  Gauge  Positions 

A  2D  side-looking  probe  was  attached  to  the  acoustic  Doppler  velocimeter 
cylinder  at  the  end  of  a  steel  rod  40  cm  in  length,  allowing  for  measurement  of 
cross-shore  and  longshore  water  velocities.  Longshore  velocities  associated  with 
three-dimensional  turbulence  were  small  compared  with  cross-shore  velocities,  which 
were  dominated  by  waves  and  wave-induced  currents.  The  3D  side-looking  probe, 
which  could  also  measure  vertical  velocities,  was  not  used  because  it  was  more  easily 
exposed  above  wave  troughs  in  shallower  water.  Thus,  as  in  Herrman  et  al.  (1997), 
only  cross-shore  velocities  were  analyzed  in  these  tests. 

Horizontal  and  vertical  locations  of  the  21  ADV  measurement  positions  are 
illustrated  in  Figure  3.16.  The  2D  probe  was  positioned  at  17  different  cross-shore 
locations  with  a  horizontal  spacing  of  0.50  m.  At  two  of  these  positions  {x  =  7.35  m 
and  X  =  9.85  m)  it  was  adjusted  vertically  to  obtain  three-point  profiles  of  the 
cross-shore  velocity  over  depth.  At  the  remaining  positions  the  probe  was  located 
as  close  as  possible  to  mid-depth.  However,  the  deep  water  probe  positions  were 
generally  above  mid-depth  to  avoid  impacts  of  wave  crests  on  the  ADV  cylinder. 
This  problem  could  have  been  avoided  by  use  of  a  longer  rod  connecting  the  probe 
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Figure  3.16:  Velocimeter  Probe  Measurement  Positions  Over  Final  Profiles 
(-  -  Test  1;  Test  2) 


and  the  ADV  cylinder,  but  the  longer  rod  might  then  have  resulted  in  excessive 
vibration  of  the  probe.  The  shallow  water  probe  positions  were  below  mid-depth  to 
avoid  exposure  of  the  probe  to  air  above  wave  trqughs. 

3.6  Repeatability  of  Tests 

After  completing  the  21  runs  for  the  free  surface  and  velocity  measurements 
for  the  two  tests,  the  repeatability  of  the  incident  and  reflected  waves  was  checked 
for  each  test.  The  separated  incident  and  reflected  wave  spectra  at  the  location  of 
wave  gauge  1  were  plotted  together  for  the  21  runs,  as  shown  in  Figure  3.17.  The 
plotted  21  spectra  for  the  incident  and  reflected  waves  in  tests  1  and  2  were  almost 
identical. 

Comparisons  were  also  made  of  the  spectral  parameters,  such  as  the  incident 
wave  spectral  peak  period,  Tp ,  the  spectral  estimate  of  the  incident  significant  wave 
height,  Hmo  =  4^/mo,  with  mo  =  zero-moment  of  the  incident  wave  spectrum,  and 
the  average  reflection  coefficient,  R  =  y(mo)r/uro,  with  (mo)r  =  zero-moment  of  the 
reflected  wave  spectrum.  In  addition,  the  separated  incident  wave  time  series  were 
analyzed  using  a  zero-upcrossing  method  to  obtain  the  time  series  parameters,  such 
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Frequency  (Hz) 


Frequency  (Hz) 


Figure  3.17:  Smoothed  Power  Density  Spectra  of  Incident  (I)  and  Reflected  (R) 
Waves  for  Test  1  and  2 


Table  3.5:  Incident  Wave  Characteristics  for  Two  Tests 


Test 

T 

-*■  s 

(s) 

Hs 

(cm) 

Tp 

(s) 

Hmo 

(cm) 

Hrms 

(cm) 

R 

1 

2.5 

20.7 

2.8 

20.3 

14.3 

0.25 

2 

1.4 

18.5 

1.6 

18.2 

12.6 

0.19 

as  the  significant  wave  period,  T*,  and  the  significant  wave  height,  Hg.  Table  3.5 
lists  the  average  values  of  Tg,  Hg,  Tp,  Hmc  ^rms,  and  R  for  the  21  runs  for  both 
tests.  Here,  is  defined  by  (2.8).  In  Tables  3.6  and  3.7,  average  wave  parameters 
are  shown  for  each  of  the  21  runs  in  tests  1  and  2,  respectively. 

The  deviations  from  the  average  values  of  Tg,  Hg,  Tp,  Hrms,  and  Hmo  in  the  21 
runs  were  generally  less  than  1%  and  did  not  exceed  4%.  The  reflection  coefficient, 
R,  varied  less  than  for  Herrman  et  al.  (1997),  with  a  nearly  constant  value  of  R  = 
0.25  for  test  1  and  a  range  of  i?  =  0.18-0.19  for  test  2.  The  waves  in  test  1  plunged 
intensely  in  a  small  region  slightly  shoreward  of  the  terrace  edge,  creating  a  large 
amount  of  suspended  sediment.  Spilling  breakers  were  more  common  in  test  2,  and 
breaking  occurred  more  continuously  over  the  bar  and  in  shallower  water  up  to  the 
swash  zone.  Visual  estimates  of  the  extent  of  breaker  and  swash  zones  are  shown 
over  the  final  equilibrium  profiles  in  Figures  3.18  and  3.19. 

The  values  for  wave  setup,  rj,  root-mean-square  wave  height,  ifrmsj  skewness, 
and  kurtosis  were  also  compared  for  gauges  1-3  in  each  run  of  both  tests.  These 
values,  shown  in  Tables  3.8  —  3.11,  are  consistent  through  each  test  and  generally 
within  small  percentage  points  of  error.  It  may  thus  be  concluded  that  the  runs  of 
tests  1  and  2  were  acceptably  repeatable. 


Table  3.6:  Time  Series  and  Spectral  Parameters  for  Test  1 


Run 

Time  Series 

Spectral 

T,(s) 

Hs(cm) 

Hmo  (cTTlJ 

Hrms  ( CTTlJ 

R 

na 

■tliM 

20.6 

mam 

14.3 

0.25 

wm 

■iliM 

20.8 

mm 

14.3 

0.25 

20.8 

20.2 

14.3 

0.25 

mm 

20.2 

■iliSl 

msm 

mm 

20.3 

liM 

msm 

20.6 

msm 

14.3 

21.0 

mm 

14.3 

jnesi 

■B 

2.5 

20.6 

14.3 

BB 

^^1 

mm 

1-10 

MM 

mm 

14.3 

0.25 

IBO 

2.5 

20.2 

14.3 

InEBl 

ina 

20.4 

14.4 

iilM 

1-13 

2.6 

20.8 

20.4 

14.4 

0.25 

1-14 

2.6 

ms^M 

20.4 

14.4 

0.25 

1-15 

2.5 

mm 

20.4 

14.4 

0.25 

1-16 

2.5 

20.6 

14.4 

1-17 

2.5 

20.8 

14.4 

iiM 

1-18 

2.6 

20.8 

14.4 

1-19 

2.5 

21.1 

20.5 

14.4 

iiMI 

1-20 

2.6 

21.1 

20.5 

14.4 

1-21 

20.6 

20.4 

14.4 

0.25 

Mean 

2.5 

20.7 

20.3 

14.3 

0.25 

Max 

2.6 

21.1 

20.5 

14.4 

0.25 

Min 

2.5 

20.2 

20.1 

14.3 

0.25 

Table  3.7:  Time  Series  and  Spectral  Parameters  for  Test  2 


Spectral 

Hs  ( cm) 

Hmo  ((^'^) 

Hrms  (cTTlJ 

R 

2-1 

■BBI 

18.2 

12.6 

Q|g| 

2-2 

DBI 

18.2 

12.6 

QQI 

2-3 

1.4 

18.0 

17.6 

12.3 

0.19 

IBB 

1.4 

18.1 

12.6 

0.19 

1.4 

18.1 

12.5 

0.19 

BB 

msM 

18.1 

IH^BB 

Qjgll 

BB 

■Em 

18.1 

QEll 

2-8 

1.4 

18.3 

18.0 

12.5 

Qg|{ 

IBB 

BB 

18.0 

QQI 

BB 

18.0 

0.19 

MW 

1.4 

18.7 

12.8 

QQ 

MW 

18.8 

18.5 

12.8 

QB 

2-13 

1.4 

18.7 

18.5 

12.8 

EDI 

2-14 

1.4 

18.8 

18.4 

12.8 

QSI 

1.4 

18.7 

18.5 

12.8 

0.19 

1.4 

18.4 

12.8 

0.19 

msm 

1.4 

12.8 

0.18 

2-18 

1.4 

18.6 

18.4 

12.8 

BEI 

1.4 

18.6 

18.3 

12.7 

ESB 

2-20 

1.4 

18.7 

18.3 

12.7 

0.18 

2-21 

1.4 

IHBi 

18.4 

12.7 

0.19 

Mean 

1.4 

18.2 

12.6 

0.19 

Max 

1.4 

18.9 

18.5 

12.8 

0.19 

Min 

1.4 

18.0 

17.6 

12.3 

0.18 

Cross-Shore  Distance  (m) 


Table  3.8:  rj  and  Hrms  for  Test  1 


Setup,  rj(cm) 

■Hrms 

Gauge  2 

Gauge  3 

HBH 

-0.16 

19^1 

MElbMI 

14.14 

14.19 

1-2 

0.11 

-0.06 

-0.06 

14.51 

14.18 

1-3 

urn 

-0.15 

■■^9! 

19991 

9^ 

■EBM 

1-4 

msm 

-0.05 

I99]9i 

mmm 

MEiM 

1-5 

IBM 

-0.06 

-0.02 

999' 

■IBM 

999 

IBM 

-0.46 

-0.36 

999 

9991 

1-7 

-0.17 

-0.24 

-0.29 

14.52 

mEmum 

mmmu 

1-8 

-0.04 

■Rl^l 

■IBi 

1999 

mtwMm 

999 

■sa 

-0.06 

f 

9]1^B 

14.48 

■EWM 

■BEBI 

1-10 

-0.16 

-0.14 

14.52 

■OBI 

9991 

1-11 

0.04 

0.08 

^■09 

mEWJfim 

999 

1-12 

0.06 

-0.12 

-0.25 

14.55 

memsm 

999 

1-13 

929 

14.66 

aiaa 

999 

1-14 

■sUfl 

99B 

■E— 

999 

■BQI 

-0.05 

■■99 

14.25 

999 

■BQ 

IBM 

14.23 

999 

1-17 

-0.29 

-0.37 

mwm 

14.25 

999 

■19 

-0.10 

9E^9 

14.23 

9E9I 

mmm 

■B1I9 

-0.55 

-0.58 

991 

14.25 

14.31 

area 

-0.03 

-0.09 

14.69 

999 

msami 

■KIM 

-0.15 

-0.33 

14.69 

999 

ibbhi 

Mean 

-0.12 

-0.17 

-0.24 

14.59 

14.20 

14.24 

Table  3.9:  Skewness  and  Kurtosis  for  Test  1 


Skewness,  s 

Kurtosis,  K 

Gauge  2 

Gauge  3 

Gauge  3 

0.66 

0.62 

0.62 

MEM 

3.67 

■Ei 

0.68 

0.63 

0.62 

3.65 

0.68 

■EH 

3.62 

1-4 

0.69 

■UH 

3.66 

1-5 

0.67 

0.64 

mmm 

1-6 

0.63 

IBH' 

masM 

0.63 

0.62 

3.92 

3.73 

3.65 

0.67 

0.64 

3.89 

3.78 

3.70 

mam 

0.66 

0.62 

3.90 

3.74 

3.66 

1-10 

0.67 

0.62 

0.61 

3.91 

3.71 

3.64 

mam 

0.64 

3.91 

3.79 

3.72 

m^m 

IBEB 

0.62 

3.86 

3.74  1 

3.64 

1-13 

0.66 

0.59 

3.86 

3.71 

3.57 

1-14 

0.66 

0.60 

3.88 

3.76 

3.63 

msm 

0.65 

0.63 

0.60 

3.87 

3.75 

1-16 

0.62 

0.60 

3.86 

3.72 

■^Hll 

1-17 

■iIM 

0.63 

0.62 

3.92 

3.75 

3.65 

1-18 

0.66 

0.63 

0.59 

3.90 

3.73 

3.59 

1-19 

0.65 

0.62 

3.93 

3.74 

3.62 

1-20 

0.66 

0.62 

0.60 

3.94 

3.59 

1-21 

0.61 

IB&ai 

liHB 

3.59 

Mean 

0.63 

0.61 

IHIH 

3.64 

Table  3.10:  rj  and  Hrms  for  Test  2 


Run 

S 

Gauge  1 

etup,  r](cm 
Gauge  2 

Gauge  3 

Gauge  1 

Hrms  (cm) 
Gauge  2 

Gauge  3 

2-1 

-0.23 

-0.22 

-0.30 

12.78 

12.62 

12.50 

2-2 

-0.08 

-0.18 

-0.12 

12.77 

12.57 

12.49 

-0.15 

-0.11 

■im 

WblEWI 

12.23 

KB 

-0.18 

-0.10 

HJgH 

12.53 

■fclkW 

-0.07 

WMW 

12.50 

■rag 

|Hinw 

-0.08 

HI 

■liifiM 

12.53 

mm 

2-7 

-0.12 

-0.05 

-0.12 

12.71 

12.51 

12.40 

\W^M 

-0.13 

-0.05 

■sail 

■OH 

mvMsm 

-0.14 

-0.11 

■lai 

■^1 

■biifcW 

m^H 

llliB 

mmmsm 

-0.06 

HSHi 

12.49 

HHH 

IBslIM 

-0.25 

BEBH 

hhi 

HHH 

2-12 

-0.17 

-0.18 

-0.24 

12.98 

Mha 

-0.13 

HEI9 

-0.19 

■feWiliW 

WfelrirW 

h&h 

mnam 

-0.14 

-0.19 

WMbVW 

12.67 

HHH 

MW 

^BMW 

-0.11 

HS&H 

MW 

MIM 

-0.15 

■bl;gw 

HHH 

2-17 

-0.19 

-0.14 

■ilH 

12.64 

MW 

-0.18 

-0.03 

HiOH 

■MfeKW 

■MTaW 

■feiriiW 

BEI 

-0.17 

-0.08 

mmm 

Mdil 

moH 

-0.19 

■rara 

MW 

IHl^l 

-0.19 

bhh 

IWfeligW 

Mean 

-0.15 

-0.11 

IHsIH 

12.59 

12.55 

Table  3.11:  Skewness  and  Kurtosis  for  Test  2 


Run 

Skewness,  s 

Kurtosis,  K 

Gauge  1 

Gauge  2 

Gauge  3 

Gauge  1 

Gauge  2 

Gauge  3 

lEH 

3.06 

3.01 

\mm 

■mfl 

3.09 

3.03 

2-3 

0.36 

3.26 

3.16 

3.12 

2-4 

0.36 

0.36 

3.20 

3.10 
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Chapter  4 


FREE  SURFACE  ELEVATIONS 

In  this  chapter,  cross-shore  surface  elevation  statistics  will  be  summarized 
and  used  to  evaluate  the  applicability  of  the  exponential  gamma  distribution  in  the 
surf  and  swash  zones.  The  wave  setup,  rj,  the  standard  deviation,  cr,  the  skewness, 
s,  and  the  kurtosis,  K  for  each  of  the  measured  time  series  of  the  free  surface 
elevation,  rj,  are  computed  under  the  assumption  of  equivalency  of  probabilistic  and 
time  averaging,  as  wajs  done  by  Herrman  et  al.  (1997).  These  325  s  time  series, 
sampled  at  the  rate  of  20  Hz,  include  the  effects  of  reflected  waves.  The  root-mean- 
square  wave  height,  i/rms,  is  determined  using  equation  (2.8).  The  normalized  free 
surface  elevation,  »/*  =  (?/  —  v)  l^i  is  then  calculated  using  the  measured  values  of 
77  and  cr.  For  each  of  the  time  series  of  ?;*,  the  probability  density  function,  /(»7*), 
is  computed  using  20  bins  in  the  range  determined  by  the  measured  maximum  and 
minimum  values  of  ?/*. 

4.1  Probability  Density  Functions  for  Gauges  1-3  for  Tests  1  and  2 

Figure  4.1  shows  the  measured  probability  density  functions  at  wave  gauges 
1-3  for  the  21  runs  in  the  two  tests.  The  exponential  gamma  distribution,  /(t;*), 
given  by  equation  (2.1)  is  also  plotted  in  Figure  4.1  for  the  measured  minimum  and 
maximum  values  of  s  for  each  test.  The  ordinate  for  the  figures  plotted  for  /(j?*)  is 
the  normalized  free  surface  elevation,  ??*,  varying  vertically.  Figure  4.1  may  be  used 
as  an  indicator  of  the  high  degree  of  repeatability  for  the  21  runs  in  each  test,  as 
the  different  runs  are  almost  indistinguishable.  The  exponential  gamma  distribution 
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Figure  4.1:  Measured  and  Computed  Probability  Distributions  for  21  Runs  at 
Wave  Gauges  1-3  for  Tests  1  and  2 


42 


changes  little  for  the  measured  range  of  s  for  each  test  and  is  in  good  agreement 
with  the  measured  distributions,  except  for  the  scatter  of  data  points  at  the  peak 
of  /(t/*).  Gauges  1-3  are  located  in  relatively  deep  water  {d  =  60  cm),  and  the 
distributions  are  approximately  Gaussian,  as  expected. 

4.2  Spatial  Variations  in  Statistical  Values  of  Free  Surface  Elevations 

The  measured  cross-shore  variations  of  fj,  Hrms,  Hnaslh,  and  s  for  tests  1 
and  2  are  shown  in  Figures  4.2  and  4.3.  The  cross-shore  coordinate  origin,  a;  =  0, 
is  located  at  gauge  1.  The  mean  water  depth  h  is  given  hy  h  =  {rj  —  Z},)  where  Zb 
is  the  bottom  elevation.  The  values  of  rf,  Hrms,  H„nslh  and  s  at  the  locations  of 
wave  gauges  1-3  are  plotted  for  all  21  runs  in  each  test.  The  swash  zones  of  the 
final  equilibrium  profiles  are  shown  as  solid  lines  in  the  upper  two  plots  of  each 
figure.  As  explained  in  Section  3.5.2,  all  the  runs  for  gauge  7  and  the  first  twelve 
runs  for  gauge  8  in  test  1  were  discarded  due  to  calibration  uncertainties.  Also, 
values  of  Hrms  and  Hrmsih  have  been  removed  for  the  first  twelve  runs  for  gauge  4  in 
test  1  because  the  gauge’s  vertical  position  shifted  after  calibration.  The  data  from 
the  repeated  measurements  at  the  locations  of  the  buried  gauges  6-10  are  plotted 
individually  without  averaging  in  order  to  indicate  that  the  visible  scatter  of  data 
points  is  largely  limited  to  the  swash  zone,  as  was  also  seen  in  Herrman  et  al.  (1997). 

As  expected,  the  wave  setup,  gradually  increases  as  waves  move  shoreward 
in  shallower  water,  becoming  tangential  to  the  beach  slope  in  the  swash  zone.  Some 
variation  is  visible  in  swash  zone  measurements  of  fj.  It  will  be  noted  that  the  most 
shoreward  measurement  of  rj  on  the  top  plot  of  Figure  4.3  appears  to  be  below 
the  average  equilibrium  sand  level  explained  in  Section  3.4.  This  is  likely  due  to 
experimental  limitations  on  the  precise  measurement  of  the  sand  bottom  position 
in  shallow  water,  as  explained  in  Section  3.5.1.  The  actual  bottom  elevation  in  the 
swash  zone  during  each  test  varied  less  than  1  cm  as  listed  in  Tables  3.3  and  3.4. 
The  statistical  analysis  also  becomes  less  reliable  in  regions  of  extremely  small  depth 
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affected  by  infrequent  wave  uprush  and  downrush.  Because  of  these  uncertainties, 
the  present  data  analysis  is  limited  to  the  region  of  mean  water  depth  h  >  1.0  cm. 
In  the  region  of  shallower  depth,  h  <  1.0  cm,  it  may  be  necessary  to  account  for 
the  probability  of  occurrence  of  dry  periods  of  zero  depth,  but  it  will  be  extremely 
difficult  to  measure  very  small  depths  above  a  movable  bottom. 

The  measured  values  of  Hrms  are  approximately  constant  in  the  region  sea¬ 
ward  of  the  bar /terrace.  They  reach  a  slight  peak  after  passing  onto  the  terrace, 
then  decrease  steadily  in  the  surf  zone  and  more  rapidly  in  the  swash  zone.  The  sand 
surface  is  again  represented  by  a  solid  line  on  these  plots,  to  illustrate  the  abrupt 
decline  in  Hrms  that  begins  near  the  still  water  shoreline.  There  was  relatively  in¬ 
tense  breaking  at  the  edge  of  the  terrace  in  test  1,  though  it  is  not  very  apparent 
from  the  relatively  constant  values  of  Hrms  seaward  of  the  surf  zone.  The  slightly 
smaller  and  shorter  waves  of  test  2  actually  show  a  decline  of  Hrms  as  they  approach 
the  breaker  zone.  The  incident  wave  conditions  have  been  listed  in  Table  3.5. 

The  measured  values  of  H„asfh  in  Figures  4.2  and  4.3  reflect  the  beach  profile 
fairly  strongly  in  deeper  water,  then  begin  to  rise  rapidly  as  they  approach  the  still 
water  shoreline.  The  landward  increase  of  Hrms/ h  is  mostly  caused  by  the  landward 
decrease  of  h.  The  decrease  of  Hrms  after  the  bar  is  rapid  enough  to  produce  a 
slight  decrease  in  Hrms/h  in  the  breaker  region  for  both  tests.  The  measured  values 
of  Hrms/h  in  the  inner  surf  zone  are  more  constant  than  the  1:16  slope  data  of 
Herrman  et  al.  (1997)  and  more  similar  to  the  field  data  on  gently  sloping  concave 
beaches  obtained  by  Thornton  and  Guza  (1982,  1983).  The  measured  values  of 
Hrms/h  approach  a  value  close  to  the  intuitive  upper  limit  of  Hrms/h  =  \/8  in  the 
swash  zone,  as  predicted  by  equation  (2.7)  with  (2.8).  The  exact  upper  limit  is 
again  difficult  to  pinpoint  because  of  the  rapid  increase  and  scatter  of  Hrms/ h  as  h 
approaches  zero,  although  the  scatter  shown  in  these  two  tests  is  less  than  that  seen 
by  Herrman  et  al.  (1997).  These  results  imply  that  the  landward  decrease  of  Hrms 


is  more  gradual  than  that  of  h  in  the  swash  zone. 

Figure  4.3  shows  that  the  skewness,  s,  increased  landward  and  then  decreased 
toward  the  still  water  shoreline,  in  accordance  with  available  data  [e.g.,  Mase  and 
Kobayashi  (1991);  Raubenheimer  et  al.  (1995)].  However,  skewness  results  in  test 
1  appear  to  have  been  significantly  influenced  by  the  intense  breaking,  which  may 
have  produced  the  dip  in  s  recorded  near  the  edge  of  the  terrace.  In  comparison 
with  the  1:16  slope  data  of  Herrman  tt  al.  (1997),  the  cross-shore  variation  of  s 
over  the  barred  and  terraced  beaches  is  less  pronounced,  indicating  the  effect  of 
the  beach  profile  on  the  irregular  wave  transformation.  In  both  tests  the  skewness 
decreases  somewhat  a.s  waves  approach  the  swash  zone,  then  increases  again  in  the 
swash  and  approaches  an  upper  limit  slightly  above  the  value  s  =  2  predicted  by 
the  exponential  distribution  given  in  equation  (2.6).  The  upper  limit  of  s  =  2 
seems  to  be  accurate  to  within  about  20%,  however.  This  trend  in  s  confirms  that 
the  probability  distribution  of  the  free  surface  elevation  is  strongly  affected  by  the 
bottom  elevation,  as  was  the  case  with  the  1:16  slope  data  of  Herrman  et  al.  (1997). 
The  free  surface  statistical  data  shown  in  Figures  4.2  and  4.3  are  listed  in  Tables  4.1 
and  4.2  in  Section  4.6. 

4.3  Measured  and  Theoretical  Probability  Density  Functions 

For  each  of  the  measurement  locations,  the  measured  probability  density 
function  of  was  compared  with  the  exponential  gamma  distribution  given  by 
equation  (2.1)  together  with  the  measured  skewness,  s.  Figures  4.4  and  4.5  show 
the  degree  of  overall  agreement  for  tests  1  and  2,  respectively.  The  comparisons 
are  separated  into  four  zones  for  each  test  on  the  basis  of  the  wave  transformation 
processes  indicated  in  Figures  3.18  and  3.19  as  well  as  the  cross-shore  variations  of 
s  shown  in  Figures  4.2  and  4.3.  The  four  zones  correspond  roughly  to  the  “shoaling 
zone”  with  little  wave  breaking,  the  “surf  zone”  with  intense  wave  breaking,  the 
“transition  zone”  between  surf  and  swash,  and  the  “swash  zone”  above  the  still 
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Figure  4.4:  Measured  and  Computed  Probability  Distributions  in  Four  Zones  for 
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Figure  4.5:  Measured  and  Computed  Probability  Distributions  in  Four  Zones  for 


water  line.  The  agreement  was  good  in  the  shoaling  zone,  and  fair  in  the  surf, 
transition,  and  swash  zones. 

The  measured  minimum  and  maximum  values  of  s  in  each  zone  are  indicated 
in  these  figures.  The  exponential  gamma  distributions  based  on  these  minima  and 
maxima  are  plotted  for  each  zone  so  that  the  data  points  would  fall  between  the 
two  theoretical  curves  if  the  agreement  were  perfect.  In  the  shoaling  and  surf  zones, 
experimental  and  theoretical  curves  maintain  a  roughly  Gaussian  shape,  except 
that  the  distribution  here  is  slightly  skewed  and  more  peaked.  For  the  Gaussian 
distribution  /(?/*)  =  0.004  at  =  ±3  and  /(»/*)  =  0.4  at  ?/*  =  0  using  equation 
(2.5).  An  anomalous  data  spike  appears  in  the  transition  zone  of  Figure  4.4;  this 
is  likely  due  to  uncertainties  involved  in  the  measurement  of  the  small  water  depth 
and  sand  bottom,  which  will  be  discussed  in  more  detail  in  the  following  section. 
In  the  swash  zone,  the  exponential  gamma  distribution  becomes  exponential  with 
s  =  2  given  by  equation  (2.6).  The  agreement  is  better  for  the  large  values  of 
7j*,  which  correspond  to  the  large  free  surface  elevations.  This  is  partially  because 
the  theoretical  distributions  for  the  large  values  of  are  not  very  sensitive  to  the 
skewness,  s. 

Separate  figures  for  each  wave  gauge  position  are  plotted  in  Figures  4.13  - 
4.34  in  Section  4.7  to  show  that  Figures  4.4  and  4.5  indeed  represent  the  overall 
agreement  for  tests  1  and  2,  respectively.  The  overall  agreement  is  very  similar  to 
the  1:16  slope  data  of  Herrman  et  al.  (1997). 

Finally,  the  measured  values  of  the  kurtosis,  K,  are  compared  with  the  the¬ 
oretical  values  of  K  based  on  equation  (2.4).  As  discussed  in  Section  2.1,  equation 
(2.4)  together  with  equation  (2.3)  may  be  used  to  derive  a  theoretical  relationship 
between  K  and  s.  This  relationship  has  been  illustrated  in  Figure  2.1  and  is  now 
redrawn  as  the  solid  line  in  Figure  4.6.  This  figure  shows  all  measured  values  of 
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(JC  —  3)  as  a  function  of  s  in  comparison  to  the  theoretical  curve,  where  it  is  as¬ 
sumed  here  that  (if  —  3)  =  6  for  s  >  2.  As  was  also  found  by  Herrman  et  al.  (1997), 
the  relation  between  K  and  s  based  on  the  exponential  gamma  distribution  tends  to 
overpredict  K  somewhat  for  s  <  2  and  can  not  explain  the  increase  of  K  for  s  >  2. 
This  suggests  that  there  may  be  a  need  for  an  empirical  fit  to  data  to  describe  this 
relationship  for  practical  applications. 

4.4  Measurements  and  Uncertainties  of  Swash  Zone  Bottom  Location 

To  determine  the  elevation  of  the  sand  surface  at  each  swash  gauge  as  accu¬ 
rately  as  possible,  the  bottom  elevation  at  each  gauge  was  measured  with  a  vernier 
scaled  pointer  gauge  for  all  runs  of  each  test.  Bottom  elevations  at  three  longshore 
positions  corresponding  to  the  profile  transect  locations  discussed  in  Section  3.4  were 
recorded  for  gauges  7-10  following  each  run.  These  were  then  averaged  to  obtain  a 
single  value  for  each  run  at  each  gauge,  which  was  used  in  calculations  of  the  mean 
water  depth,  /i,  and  other  statistical  parameters.  This  averaging  is  consistent  with 
the  average  equilibrium  profile  discussed  in  Section  3.4. 

Although  there  was  significant  motion  of  the  sand  particles  during  each  400  s 
wave  run,  the  surface  measured  before  and  after  runs  varied  little.  The  variations  in 
the  bottom  depth  or  elevation  at  gauges  7-10  over  tests  1  and  2,  tabulated  earlier  in 
Tables  3.3  and  3.4,  are  now  shown  in  Figure  4.7.  Average  bottom  elevation  at  the 
four  gauges  varied  by  less  than  1  cm  over  the  entire  course  of  test  1 ,  excluding  the 
repositioning  of  these  gauges  between  runs  12  and  13.  This  variation  was  greater  at 
gauge  7  during  test  2,  which  initially  was  located  slightly  below  the  swash  zone  and 
on  the  edge  of  a  region  of  ripples.  Except  for  this  irregularity,  the  average  bottom 
elevation  measured  at  each  swash  gauge  over  test  2  also  varied  by  less  than  1  cm, 
excluding  the  repositioning  of  the  gauges  between  runs  10  and  11.  The  longshore 
deviations  from  these  averages  were  less  than  0.6  cm. 
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Figure  4.7:  Measured  Depth/Elevation  of  Sand  Bottom  at  Gauges  7—10  During 
Tests  1  and  2 
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Figure  4.8:  (A)  Uncorrected  Time  Series  of  Free  Surface  at  Swash  Gauge  9  for 
Run  5  of  Test  2;  (B)  Corrected  Time  Series 


The  bottom  elevation,  zi,  in  the  swash  zone  directly  affects  the  instantaneous 
water  depth  h  =  {i]  -  Zi),  where  the  free  surface  elevation  rj  was  measured  by  the 
swash  gauges.  In  the  swash  zone,  h  can  be  very  small  and  the  small  variation  of 
zi  influences  the  value  of  h  significantly.  When  the  sand  surface  was  exposed  at  a 
given  swash  gauge,  the  gauge  reading  would  correspond  to  the  bottom  elevation, 
zi),  only  if  the  sand  were  completely  saturated.  Otherwise  the  reading  would  drop 
below  this  value  to  reflect  the  lower  water  level  beneath  the  sand  surface.  Because 
of  uncertainty  about  how  the  gauge  capacitance  would  be  affected  by  partially  sat¬ 
urated  sand,  it  was  decided  not  to  attempt  to  measure  the  water  level  below  the 
sand  surface.  In  the  time  series  of  surface  elevations  in  the  swash  zone,  all  gauge 
readings  below  the  measured  sand  bottom  value  were  set  equal  to  that  value.  This 
“cutoff”  procedure  is  shown  for  a  typical  swash  gauge  time  series  in  Figure  4.8. 

In  Figure  4.8(A)  the  measured  sand  bottom,  Zh,  is  indicated  by  the  horizontal 
line  drawn  at  0.004  m.  The  gauge  readings  in  this  400  s  time  series  fall  below  the 
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sand  bottom  many  times,  occasionally  by  distances  approaching  2  cm.  Analysis  of 
time  series  such  as  this  required  that  all  values  below  this  line  be  reassigned  the 
measured  value  of  Zb,  as  shown  in  Figure  4.8(B).  The  adjusted  time  series  of  this 
form  were  then  analyzed  to  obtain  the  swash  zone  statistical  values  presented  in  this 
chapter. 

Several  factors  may  have  contributed  to  the  difference  between  measured 
bottom  position  and  minimum  swash  gauge  readings  illustrated  in  Figure  4.8.  The 
suspension  of  large  quantities  of  sediment  during  wave  uprush  and  downrush  in  the 
swash  zone  lowers  the  apparent  bottom  elevation.  The  water  was  consistently  very 
cloudy  during  the  runs,  and  it  was  not  possible  to  measure  sand  elevations  under 
flowing  water.  Since  the  swash  gauges  were  calibrated  under  quiescent  water,  it 
is  not  certain  whether  gauge  readings  exactly  corresponded  to  the  water  elevation. 
Additionally,  because  of  some  longshore  variation  of  the  beach  profile  the  three-point 
averaged  Zb  values  did  not  always  reflect  the  exact  bottom  elevation  at  the  location 
of  each  gauge.  However,  this  difference  was  generally  less  than  0.3  cm.  In  light  of 
these  uncertainties,  the  measurement  of  water  depths,  h,  less  than  about  1  cm  may 
not  have  been  very  accurate. 

The  probability  distributions  shown  in  Figure  4.4  included  an  anomalous 
group  of  data  points  in  the  third  plot,  corresponding  to  approximately  ?/*  =  -1.5. 
These  values,  which  can  also  be  seen  clearly  in  Figure  4.22,  imply  that  the  very 
small  water  depth  occurred  more  often  than  predicted  by  the  exponential  gamma 
distribution.  Lowering  the  specified  bottom  elevation,  Zb,  would  increase  the  water 
depth,  h,  and  could  eliminate  the  anomalous  data  points.  This  change  is  not  made 
herein,  in  order  to  clearly  indicate  the  uncertainties  and  difficulties  in  measuring  the 
small  water  depth  above  the  movable  bed.  Similar  spikes  also  appear  in  the  fourth 
plot  of  Figures  4.4  and  4.5,  but  these  appear  to  be  an  extension  of  the  exponential 
distribution  for  s  =  2.  The  uncertainties  in  bottom  elevation  and  the  cutoff  line 
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discussed  above  directly  affect  the  position  of  these  spikes  in  each  figure.  The  value 
rj^  =  —1  for  the  exponential  distribution  corresponds  to  7]  =  Zb,  as  explained  in 
relation  to  equation  (2.6).  Consequently,  the  measured  probability  distributions  in 
the  vicinity  of  ?/*  =  —  1  for  the  swash  zone  may  not  be  very  accurate. 


4.5  Wave  Spectral  Transformation  in  Cross-Shore  Direction 

A  preliminary  examination  was  conducted  into  the  transformation  of  the  total 
(incident  -|-  reflected)  wave  spectra  recorded  at  different  cross-shore  locations.  The 
object  of  this  examination  was  to  detect  evidence  of  the  existence  of  standing  waves 
over  the  terrace  region  in  test  1  or  test  2.  The  first  mode  of  such  a  standing  wave, 
illustrated  for  test  1  in  the  bottom  diagram  of  Figure  4.9,  would  have  a  node  located 
over  the  edge  of  the  terrace  with  antinodes  located  both  offshore  and  near  the  still 
water  shoreline.  The  wavelength  of  such  a  standing  wave  would  be  approximately 
L  ~  4/,  where  I  ~  5.66  m  is  the  distance  from  t.he  terrace  edge  to  the  still  water 
shoreline  in  test  1.  Thus  for  test  1,  L\  —  22.64  m.  For  test  2,  with  I  —  6.91  m,  we 
find  Li  ~  27.64  m. 

Figures  4.9  and  4.10  show  wave  spectra  measured  in  the  vicinity  of:  (1)  the 
offshore  antinode  position;  (2)  the  node  position;  (3)  a  mid-terrace  location;  and  (4) 
the  shoreline  antinode  position  for  a  first  mode  standing  wave  in  test  1  and  test  2, 
respectively.  Were  such  a  standing  wave  present,  it  would  be  identifiable  as  a  peak 
in  spectra  for  positions  (1)  and  (4)  and  a  trough  in  the  spectrum  at  position  (2), 
all  occurring  at  the  same  low  frequency.  Such  consistent  low  frequency  peaks  and 
troughs  are  not  apparent  in  either  of  these  two  figures. 

Mode  frequencies  are  estimated  by  the  simple  relationship 


C  _  \[^t 

Ti  ~  4//(2f  -  1) 


for 


1,2,... 


(4.1) 


where  ht  =  mean  depth  on  the  terrace  ~  23  cm,  and  linear  long  wave  theory  is 
applied  since  Li  >  'Zbhi  ~  6  m  when  i  =  1  or  2.  From  equation  (4.1),  we  find  that 
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Position  1  {2.60m) 


Position  2  (7.85m) 


Figure  4.9:  Unsmoothed  Wave  Spectra  at  Four  Cross- Shore  Locations  Near  Nodes 
and  Antinodes  of  First  Mode  Standing  Wave,  Test  1 


Figure  4.10:  Unsmoothed  Wave  Spectra  at  Four  Cross-Shore  Locations  Near 
Nodes  and  Antinodes  of  First  Mode  Standing  Wave,  Test  2 


for  test  1,  /i  =  0.07  Hz  and  /2  =  0.20  Hz,  while  for  test  2,  /i  =  0.05  Hz  and  /2  = 
0.16  Hz. 

In  the  top  half  of  Figures  4.11  and  4.12  the  wave  spectra  at  position  (4) 
from  tests  1  and  2  are  expanded  to  examine  more  closely  the  locations  of  spectral 
peaks  at  the  respective  shorelines.  Two  vertical  dashed  lines  have  been  added  to 
each  spectral  plot  in  order  to  indicate  the  standing  wave  frequencies  for  the  first 
(/i)  and  second  (/2)  modes.  The  envelopes  of  both  modes  are  displayed  with  the 
corresponding  equilibrium  beach  profile  in  the  lower  half  of  each  figure.  Although 
there  is  a  localized  peak  around  the  second  mode  frequency  in  Figures  4.11  and  4.12, 
these  standing  waves  are  not  dominant.  The  lack  of  clear  evidence  of  standing  waves 
of  the  first  or  second  mode  in  either  test  indicates  that  the  barred  and  terraced 
profiles  in  test  1  and  2  are  not  caused  by  standing  waves. 
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Depth  (m)  cm*cm*s 


Position  4  (13.51nn) 


Figure  4.11:  Expanded  Wave  Spectrum  at  Still  Water  Shoreline,  Test  1  (Standing 
Wave  First  Mode  Frequency  /i  ^  0.07  Hz,  Second  Mode  Frequency 
/2  ~  0.20  Hz) 


Depth  (m)  cm*cm’"s 


Position  4  (13.51m) 


Frequency  (Hz) 
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Figure  4.12:  Expanded  Wave  Spectrum  at  Still  Water  Shoreline,  Test  2  (Standing 
Wave  First  Mode  Frequency  fi  ~  0.05  Hz,  Second  Mode  Frequency 
/2  ~  0.16  Hz) 


4.6  Complete  Free  Surface  Statistical  Data  For  Tests  1  and  2 

Complete  free  surface  statistical  data  is  given  below  for  each  gauge  location 
for  tests  1  and  2.  Each  gauge  is  indexed  by  its  x  position  and  the  corresponding 
water  depth,  where  Zf,  =  —d.  For  gauges  1,  2,  and  3  located  at  a:  =  0,  0.38,  and 
0.50  m,  respectively,  the  measured  values  of  IJ,  ifnns,  -s,  and  K  for  the  21  runs  have 
been  listed  in  Tables  3.8  -  3.11.  The  values  listed  in  Tables  4.1  and  4.2  are  the 
average  values  for  these  three  gauges  over  the  21  runs. 


Table  4.1:  Test  1  Gauge  Data 


X 

(m) 

(cm) 

V 

(cm) 

h 

(cm) 

Hrms 

(cm) 

s 

Kurtosisj  K 

data 

theory 

0.00 

-60.0 

-0.12 

59.88 

3.90 

3.87 

0.38 

-60.0 

-0.17 

59.83 

14.20 

0.63 

3.73 

3.78 

IlilECT 

-60.0 

-0.24 

59.76 

14.24 

0.61 

3.64 

3.74 

57.99 

— 

0.84 

4.31 

4.35 

2.85 

-58.1 

-0.23 

57.87 

13.36 

0.82 

4.14 

mill' 

mm 

56.94 

— 

0.80 

4.09 

mwjum 

BiliEl 

0.82 

4.14 

3.60 

-56.0 

-0.16 

55.84 

— 

0.89 

4.42 

55.39 

0.94 

4.65 

liiil 

-54.1 

-0.40 

53.70 

— 

-52.8 

-0.32 

0.92 

mmm 

mmi 

51.40 

— 

BlEa 

-50.1 

0.00 

50.10 

13.56 

-48.7 

-0.48 

48.22 

— 

5.35 

-47.2 

-0.42 

46.78 

13.37 

0.92 

4.47 

-45.1 

■nia 

44.92 

— 

-42.6 

Bina 

0.87 

4.18 

■EB 

6.10 

-40.3 

-0.33 

39.97 

— 

0.87 

4.22 

iBIB 

maem 

■iifigl 

■waii 

HSI 

IgllEIIl 

Ktmiw 

II^H 

0.92 

4.39 

4.63 

6.85 

-31.7 

-0.30 

31.40 

13.94 

0.91 

4.31 

4.58 

Table  4.1,  continued 


■ 

Zh 

V 

h 

H'prns 

Ml 

(cm) 

(cm) 

(cm) 

(cm) 

7.10 

B^l 

— 

0.77 

3.97 

K^l 

7.10 

BeWrl 

■agl 

27.28 

— 

0.78 

4.02 

4.18 

MIM 

-27.7 

-0.46 

■HI 

■riKM 

-23.6 

-0.26 

0.76 

3.88 

KSfl 

mism 

23.41 

14.30 

0.77 

3.86 

mm 

■QKfiii 

23.22 

14.27 

7.60 

ren 

20.76 

15.79 

0.66 

3.70 

3.86 

■rliiSB 

-21.3 

-0.11 

0.88 

4.04 

mm 

renal 

-22.5 

-0.87 

21.62 

15.05 

m^m 

8.35 

-23.1 

-0.56 

13.71 

0.86 

3.98 

IQB 

-23.3 

-0.21 

MtillM 

13.15 

0.96 

4.03 

■a 

B:i:re 

23.10 

12.86 

renal 

-22.2 

-0.02 

22.13 

12.49 

urn 

BIgfel 

-21.4 

0.06 

miji 

gMBl 

BBI 

9.60 

-21.1 

0.29 

21.39 

11.07 

0.91 

3.70 

retsM 

MWii 

■a 

mam 

0.97 

4.02 

4.79 

21.16 

10.87 

Ij^H 

n^ii 

■BB 

-20.7 

0.42 

■ikU 

KiiBl 

9.85 

20.93 

10.90 

0.88 

3.59 

mna 

-19.7 

0.49 

20.21 

0.91 

3.62 

18.87 

10.65 

0.86 

3.35 

4.44 

EIEEI 

16.96 

10.60 

0.95 

3.74 

mm\ 

■lil^ 

-15.5 

0.48 

15.98 

0.89 

3.64 

IBB 

11.55 

0.69 

3.19 

11.55 

mdwM 

■Elia 

0.66 

KWBM 

■tiM 

11.55 

Bfciil 

■Elrfil 

Bfl 

3.90 

11.55 

■HEl 

0.68 

■Elrf;l 

11.55 

mam 

■fclilil 

0.68 

3.11 

Idi 

-13.9 

1.23 

15.13 

mi 

3.95 

lairei 

-13.9 

0.80 

■Erei 

0.70 

3.21 

■aBEM 

imrei 

iniTOa 

|E^ 

lEI^ 

11.55 

-13.9 

0.70 

14.60 

8.81 

0.67 

3.15 

3.88 
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Table  4.1,  continued 


X 

Zh 

V 

h 

Hrms 

5 

Kurtosisj  K 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

data 

theory 

11.55 

-13.9 

0.81 

14.71 

8.81 

0.69 

3.27 

3.94 

11.55 

-13.9 

0.66 

mil 

11.55 

-13.9 

lEIilrl 

WMM 

mui 

-11.2 

1.02 

12.22 

msi 

liWl 

3.36 

ggg 

-11.2 

0.65 

11.85 

KliCTi 

BWl 

3.37 

gim 

12.25 

-11.2 

0.94 

12.14 

8.50 

0.73 

3.36 

4.05 

-11.2 

0.68 

11.88 

0.75 

3.36 

4.09 

-11.2 

0.84 

12.04 

mEi 

0.73 

3.34 

4.05 

■ilgfli 

12.10 

8.45 

0.74 

3.41 

4.08 

IlMfefci 

■fjQI 

12.00 

8.44 

0.70 

3.30 

3.98 

12.25 

-11.2 

0.84 

12.04 

8.42 

0.71 

3.31 

3.99 

,40 

0.83 

12.03 

msii 

Em 

giig 

Ifcfifci 

1.27 

3.57 

mu 

mKvm 

13.65 

0.67 

MiM 

0.83 

13.65 

1.10 

QQI 

0.70 

ItMa 

13.65 

-2.3 

0.56 

2.86 

6.67 

0.88 

3.55 

4.48 

QQI 

0.54 

2.84 

0.92 

3.55 

4.64 

Ifcfifci 

0.56 

2.96 

0.82 

3.37 

4.29 

-2.5 

■IMM 

3.23 

■HIM 

0.79 

mg 

4.22 

mail 

-2.4 

3.14 

0.76 

4.12 

13.65 

-2.6 

0.43 

3.03 

6.79 

0.81 

3.39 

4.28 

■hWil 

2.1 

4.11 

2.01 

mu 

1.40 

igi 

gim 

■fclrifel 

2.0 

3.97 

1.97 

1.40 

misi 

■dEf 

13.79 

1.9 

3.74 

6.78 

urn 

4.32 

6.37 

13.79 

1.8 

3.71 

mg 

6.75 

igi 

4.45 

6.40 

13.79 

1.7 

3.76 

2.06 

6.94 

1.34 

4.33 

6.20 

moisi 

3.44 

6.69 

igi 

misi 

6.45 

■fclriil 

mm 

3.69 

mu 

6.96 

lig 

imi 

6.21 

13.79 

1.7 

3.74 

2.04 

mu 

1.34 

4.32 

13.79 

1.6 

3.63 

2.03 

mu 

1.32 

4.25 

giiMI 

13.79 

1.5 

3.49 

1.99 

7.03 

1.34 

4.25 

6.17 

13.79 

3.40 

m&i 

imu 

1.34 

gm 

mu 

13.79 

ma 

3.33 

■llika 

mu 

1.34 

ggg 
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Table  4.1,  continued 


X 

(m) 

Zb 

(cm) 

7 

(cm) 

h 

(cm) 

Hrms 

(cm) 

5 

Kurtosis,  K 

data 

theory 

13.93 

4.0 

5.21 

1.21 

4.64 

1.88 

6.97 

8.48 

4.84 

dE 

dsll 

1.97 

7.39 

3.6 

4.99 

Idl 

Edi 

■iMl 

ddl 

13.93 

3.6 

4.88 

1.28 

4.76 

1.84 

7.08 

8.31 

■kMI 

MifEMI 

■ifiMI 

Od 

1.74 

6.66 

7.88 

mum 

da 

Ed 

MEM 

1.74 

6.24 

7.90 

3.4 

ESI 

1.51 

Edi 

1.61 

5.86 

7.31 

3.2 

HE 

1.21 

eei 

1.71 

5.99 

7.77 

13.93 

3.1 

4.46 

1.36 

4.75 

7.43 

14.07 

7.4 

8.31 

■tIilB 

dS 

2.37 

ddl 

14.07 

7.4 

8.37 

■tliM 

2.35 

8.65 

■emli 

msm 

0.92 

dS 

2.30 

8.56 

9.00 

■Elilrl 

8.20 

Idl 

9.00 

14.07 

7.1 

8.11 

1.01 

3.64 

mm 

■EBMI 

Ed 

Qg| 

2.17 

7.68 

mm 

lEiilrl 

dE 

Idl 

3.75 

9.00 

IlliU 

S3 

dE 

1.08 

3.79 

dSdl 

■Elilrl 

dE 

1.23 

3.80 

2.08 

7.18 

eeiI 

14.07 

7.0 

8.06 

1.06 

2.07 

6.96 

9.00 

dia 

■ilM 

dE 

1.14 

Ed 

dE 

9.00 

mm 

d£l 

dE 

9.00 

14.21 

10.0 

10.64 

d£l 

10.34 

9.00 

14.21 

9.9 

10.49 

■iligl 

2.71 

10.86 

9.00 

ilEldl 

9.8 

10.48 

eei 

2.71 

10.97 

dd 

■Ed 

9.8 

10.50 

Ed 

mmm 

9.00 

14.21 

9.6 

10.40 

Ed 

2.63 

10.45 

9.00 

14.21 

9.5 

■lllill 

■IMHM 

msm 

2.77 

11.30 

9.00 

14.21 

■Itfcfcl 

EEE 

9.00 

14.21 

9.3 

10.30 

1.00 

2.90 

2.64 

10.58 

9.00 

14.21 

9.2 

9.98 

0.78 

2.99 

2.52 

9.60 

9.00 
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Table  4.2:  Test  2  Gauge  Data 


X 

(m) 

Zb 

(cm) 

V 

(cm) 

h 

(cm) 

Hrms 

(cm) 

s 

Kurtosis,  K 

data 

theory 

0.00 

-60.0 

-0.15 

59.85 

12.81 

0.39 

3.18 

3.31 

0.38 

-60.0 

EDII 

BH 

0.50 

-60.0 

■iHIBi 

2.60 

-57.2 

12.53 

2.85 

-56.9 

-0.23 

56.70 

12.50 

0.23 

2.95 

Ba 

3.10 

-55.0 

-0.21 

Ba 

3.35 

QOI 

Ba 

3.60 

-52.8 

-0.13 

52.71 

12.06 

0.24 

3.00 

3.11 

3.85 

51.34 

11.96 

3.19 

-49.8 

49.75 

12.17 

3.21 

glia 

-48.0 

BQI 

0.37 

3.04 

3.28 

4.60 

-46.1 

-0.26 

Bsai 

4.85 

Biisa 

KBS 

0.43 

3.09 

5.10 

HQI 

41.23 

12.14 

0.43 

3.07 

3.36 

-0.13 

38.15 

0.47 

3.16 

3.44 

01^9 

mm 

34.50 

3.42 

5.85 

-31.3 

-0.20 

31.08 

12.07 

0.52 

3.26 

IBBI 

6.10 

-28.2 

-0.18 

0.52 

3.30 

■tiMi 

6.35 

3.70 

6.60 

0.63 

3.39 

3.78 

6.85 

-22.2 

-0.23 

■fellfcg 

iigEa 

4.10 

IBEI 

-22.3 

-0.21 

22.09 

11.69 

4.34 

7.10 

iioa 

4.37 

7.10 

-22.3 

-0.23 

22.07 

12.05 

0.82 

3.51 

7.35 

-23.6 

-0.09 

regal 

BIB 

0.95 

3.75 

7.35 

Inra 

reifcjii 

BE£| 

0.91 

3.58 

4.60 

7.35 

-23.6 

-0.18 

23.42 

11.24 

0.89 

3.65 

4.51 

Table  4.2,  continued 


X 

(m) 

(cm) 

(cm) 

h 

(cm) 

Hrms 

(cm) 

S 

Kurtosis,  K 

data 

theory 

7.60 

-24.5 

-0.05 

24.43 

10.75 

4.91 

7.85 

-24.5 

-0.01 

24.50 

10.32 

1.02 

3.72 

4.97 

8.10 

24.52 

9.88 

1.00 

3.68 

4.89 

8.35 

-24.5 

0.06 

24.57 

9.65 

0.96 

3.44 

8.60 

24.13 

9.47 

0.93 

3.41 

-23.4 

0.21 

23.59 

9.38 

4.56 

MKm 

^^1 

3^1 

0.89 

3.25 

BOB 

-21.8 

0.19 

mm 

0.88 

3.24 

BEifl 

-20.9 

0.24 

BHI 

■;«iM 

031 

3.24 

_ 

4.44 

MM 

4.49 

21.09 

8.93 

0.89 

3.29 

4.51 

-19.9 

0.18 

20.08 

8.82 

0.90 

3.37 

BBBI 

8.91 

0.89 

3.26 

9.85 

3QI 

3[g 

£31 

£31 

4.54 

10.10 

-19.2 

0.22 

■aECTi 

4.64 

ll■lll^ali 

-17.5 

0.26 

17.74 

8.52 

0.91 

3.51 

4.60 

10.85 

-16.9 

8.43 

0.96 

3.58 

3^3 

11.55 

-15.6 

1013 

7.86 

3E!l!l 

IfcfcEl 

7.88 

1.03 

3.65 

5.00 

-15.6 

0.31 

■fcKiH 

7.78 

1.03 

3.73 

5.00 

11.55 

15.85 

7.78 

133 

MM 

—ill 

3.66 

BIB 

11.55 

-15.6 

0.34 

■fclZl 

IMBSa 

inEEl 

BEBI 

11.55 

-15.6 

0.26 

15.86 

7.76 

bsb 

11.55 

-15.6 

15.83 

7.73 

1.02 

3.65 

4.97 

11.55 

■waiH 

319 

Bail 

11.55 

-15.6 

0.27 

IBM 

Table  4.2,  continued 


Zb 

(cm) 

7 

(cm) 

h 

(cm) 

H-rms 

(cm) 

5 

Kurtosis,  K 

data 

theory 

12.25 

Bmi 

14.29 

7.66 

4.82 

my 

-13.9 

0.35 

imii 

BSfl 

1.00 

3.57 

Ifelifci 

lEMlI 

mu 

0.98 

3.60 

■mi 

■Eilill 

1.00 

3.57 

■mi 

-13.9 

0.43 

14.33 

7.67 

■mi 

12.25 

lEiW 

0.96 

3.49 

4.75 

imil 

-13.9 

0.27 

iim 

0.97 

3.46 

4.79 

-13.9 

0.41 

|EM| 

0.97 

3.54 

4.78 

12.25 

-13.9 

0.33 

14.23 

7.70 

0.97 

3.55 

4.78 

1.01 

3.60 

■smi 

imi' 

0.32 

9.12 

7.90 

■■11 

iim 

-9.8 

0.53 

3.78 

-7.7 

0.48 

8.18 

0.68 

3.22 

3.92 

13.23 

W;M 

3.86 

13.23 

-9.7 

0.36 

■lilia 

mm 

■wara 

Kija 

3.88 

13.23 

-9.5 

0.35 

7.83 

0.66 

3.31 

3.87 

ms 

ms 

0.68 

3.31 

3.91 

■fcliSl 

azfia 

3.74 

13.23 

-7.9 

0.18 

8.08 

8.20 

0.62 

3.16 

13.23 

-7.7 

0.22 

ms 

8.22 

0.58 

3.14 

13.37 

-5.0 

0.49 

ms 

6.96 

0.86 

3.49 

■sm 

BS 

ms 

ms 

0.83 

3.40 

4.33 

■ilia 

MM 

Bjm 

4.40 

13.37 

-5.6 

0.41 
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4.35 

13.37 

-5.6 
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3.39 

■m 
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mm 
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ms 
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ms 

ms 
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6.06 
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im 

ms 
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mm 
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-5.8 
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Table  4.2,  continued 


H 

Zb 

(cm) 

V 

(cm) 

h 

(cm) 

Hrms 

(cm) 

5 

Kurtosisj  K 

theory 

|gg| 

-4.9 

0.65 

5.55 

4.36 

■fciitl 

-4.6 

0.79 
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3.78 
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BM\ 
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nn 

0.85 
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mm 
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Bn 
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i^n 
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nn 
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13.65 
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■W«l^ 

MSM 
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3.86 

■m 

13.65 

-3.0 

mmm 

B|gB 
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3.85 

13.65 
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1.31 

4.31 

5.95 

4.18 

— 

13.65 

BSI 

un 

Bn 

bsei 

0.73 

3.84 

4.05 
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1.22 
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■WEB 

mi 

ms 
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■fcfgKi 

-3.1 

■KOT 

Bliltl 

nn 

0.74 

3.76 

Bsn 

13.65 

-3.0 

■ICT 

4.24 
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3.86 

-3.0 

■m 

Bn 

nB 
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3.85 

llfcfdil 

Mmm 

nB 

llg^ 

m 
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3.83 
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Bn 

Bn 

^Bl 

■HkM 

IlfcyCT 

mam 

2.64 

iBa 

4.62 

1.13 

4.54 

5.36 

13.79 

0.0 

muniM 

imi 

1.08 

4.25 

nn 

13.79 

0.0 

■ill;! 

lEOI 

bee 

11^ 

0.4 

2.64 

2.24 

11^ 

0.4 

2.64 

■SB 

BlgB 

13.79 

lifckiM 

1.60 

4.39 

1.36 

5.15 

6.29 

IQQI 

0.2 

2.54 

iBn 

IggBl 

IlhlfiCT 

0.1 

2.36 

IldMiW 

1^11 

13.79 

0.0 

2.42 

2.42 

4.89 

1.17 

4.71 

5.53 
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Table  4.2,  continued 


Zb 

V 

h 

Hrms 

(cm) 

(cm) 

(cm) 

(cm) 

s  Kurtosis,  K 


data 


13.93 

1.9 

3.78 

1.88 

4.54 

1.48 

6.03 

13.93 

2.0 

2.02 

1.43 

6.12 

14.07 

4.9 

0.52 

■KM 

1.78 

6.68 

8.04 

lEliM 

■WiH 

KtIrfrK 

2.02 

5.25 

2.08 

1.80 

7.71 

14.07 

4.6 

5.40 

0.80 

2.10 

1.58 

6.05 

■EKMMPM 

5.69 

m 

.69 

2.07 

Q3 

5.96 

BBI 

■iliaggM 

5.64 

ffi 

.54 

1.98 

IQI 

7.38 

BE! 

■Htia 

5.61 

0.61 

1.70 

6.68 

mifim 

■Elia 

4.7 

5.40 

0.70 

1.68 

6.85 

14.07 

4.7 

5.44 

0.74 

2.18 

1.50 

5.55 

14.07 

4.5 

5.37 

1.58 

6.13 

IBi 

14.21 

2.09 

7.99 

9.00 

14.21 

6.78 

2.13 

8.16 

■iMIIlM 

14.21 

■aw 

6.75 

2.08 

8.08 

,4.2: 


14.21  6.2 


6.2 


6.2 


6.1 


14.21  6.1 


0.47 

1.81 

gdl:feK 

0.62 

■KfcM 

2.03 

7.98 

0.63 

■KM 

2.07 

8.12 

6.78 

2.00 

7.79 

6.68 

1.85 

2.12 

8.53 

1.92 

2.03 

7.86 

0.47 

1.85 

8.67 

4.7  Measured  and  Computed  Probability  Distributions  at  Each  Gauge 
Position 

On  the  following  pages  Figures  4.13  -  4.34  show  individual  plots  of  measured 
versus  computed  probability  distributions  of  the  free  surface  at  each  measurement 
position  for  tests  1  and  2,  excluding  the  distributions  at  gauges  1-3,  which  were 
already  shown  in  Figure  4.1.  Note  that  repeated  measurements  at  a  single  location 
for  the  gauges  buried  in  the  sand  and  for  multiple  measurements  at  velocity  profile 
locations  are  plotted  together. 
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Figure  4.13:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  2.60,  2.85,  3.10,  and  3.35  m  for  Test  1 
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Test  1 


Figure  4.15:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  4.60,  4.85,  5.10,  and  5.35  m  for  Test  1 


Test  1 


Figure  4.16:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  5.60,  5.85,  6.10,  and  6.35  m  for  Test  1 


Test  1 


Figure  4.17:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  6.60,  6.85,  7.10,  and  7.35  m  for  Test  1  (x  =  7.10  and  x  =  7.35 
are  Velocity  Profile  Locations) 


Test  1 


Figure  4.18:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  7.60,  7.85,  8.10,  and  8.35  m  for  Test  1 


Test  1 


Figure  4.19:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  8.60,  8.85,  9.10,  and  9.35  m  for  Test  1 


Test1 


Figure  4.20:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  9.60,  9.85,  10.10,  and  10.35  m  for  Test  1  (x  =  9.60  and  x  =  9.85 
are  Velocity  Profile  Locations) 
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Probability  Density  Function  f 


Figure  4.21:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  10.60,  10.85,  11.55,  and  12.25  m  for  Test  1 
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Test  1 


Figure  4.22:  Measured  and  Computed  Probability  Distributions  at  Position 
X  =  13.65  m  for  Test  1  (Positions  x  =  13.23,  13.37,  and  13.51  were 
discarded) 
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Figure  4.23:  Measured  and  Computed  Probability  Distributions  at  Positions 
x  =  13.79,  13.93,  14.07,  and  14.21  m  for  Test  1 


Test  2 


Figure  4.24:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  —  2.60,  2.85,  3.10,  and  3.35  m  for  Test  2 


Test  2 


Figure  4.25:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  3.60,  3.85,  4.10,  and  4.35  m  for  Test  2 


Test  2 


Figure  4.26:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  4.60,  4.85,  5.10,  and  5.35  m  for  Test  2 


Test  2 


Figure  4.27:  Measured  and  Computed  Probability  Distributions  at  Positions 
a:  =  5.60,  5.85,  6.10,  and  6.35  m  for  Test  2 


Test  2 


Figure  4.28:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  6.60,  6.85,  7.10,  and  7.35  m  for  Test  2  (x  =  7.10  and  x  =  7.35 
are  Velocity  Profile  Locations) 
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Test  2 


Figure  4.29:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  —  7.60,  7.85,  8.10,  and  8.35  m  for  Test  2 


Test  2 


Figure  4.30:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  8.60,  8.85,  9.10,  and  9.35  m  for  Test  2 


Test  2 


Figure  4.31:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  9.60,  9.85,  and  10.10  cm  for  Test  2  (a;  =  9.60  and  x  =  9.85 
are  Velocity  Profile  Locations;  Data  at  a;  =  10.35  was  discarded) 
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Figure  4.32:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  10.60,  10.85,  11.55,  and  12.25  m  for  Test  2  (Run  20  for  gauge  6 
at  X  =  12.25  m  was  discarded) 
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Test  2 


Figure  4.33:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  13.23,  13.37,  13.51,  and  13.65  m  for  Test  2 


Test  2 
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Figure  4.34:  Measured  and  Computed  Probability  Distributions  at  Positions 
X  =  13.79,  13.93,  14.07,  and  14.21  m  for  Test  2 


Chapter  5 


CROSS-SHORE  HORIZONTAL  VELOCITIES 

In  the  present  chapter,  cross-shore  velocity  statistics  will  be  summarized  and 
compared  to  the  predictions  of  the  exponential  gamma  distribution  in  surf  and  swash 
zones,  as  was  done  for  free  surface  elevations  in  chapter  4.  Also  included  in  this 
analysis  will  be  an  examination  of  the  variation  of  these  statistics  over  depth  and  an 
attempt  to  synchronize  free  surface  and  velocity  time  series.  Velocity  measurements 
were  made  at  17  cross-shore  locations,  with  two  3-point  vertical  profiles  taken  at 
X  =  7.35  m  and  x  =  9.85  m.  The  mean  velocity,  u,  the  standard  deviation,  cr^, 
the  skewness,  s^,  and  the  kurtosis,  K^,  are  calculated  for  each  time  series  of  the 
cross-shore  velocity,  u,  measured  at  the  approximate  mid-depth  below  SWL.  These 
are  computed  under  the  assumption  of  equivalency  of  the  probabilistic  and  time 
averaging,  in  the  same  manner  as  Herrman  et  al.  (1997).  The  horizontal  velocity, 
u,  is  taken  as  positive  landward,  so  that  the  undertow  velocity,  u,  is  negative. 
The  normalized  horizontal  velocity,  =  (u  — u)/(7u,  is  then  calculated  using  the 
measured  values  of  u  and  (7„.  The  probability  density  function,  /(«»),  for  each  time 
series  of  u*  is  determined  and  compared  with  the  exponential  gamma  distribution. 
The  theoretical  distribution  is  given  by  equations  (2.1)  -  (2.3),  where  and  s  must 
be  replaced  by  u*  and  s^,  respectively.  Skewness  and  kurtosis  values  for  u  are  the 
same  as  those  for  u*. 
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5.1  Horizontal  Velocity  and  Linear  Theory 

In  obtaining  approximate  relationships  between  the  free  surface  and  horizon¬ 
tal  velocity  statistics,  local  nonlinearity  is  neglected  and  linear  progressive  long- wave 
theory  is  assumed  to  be  approximately  valid  locally  even  inside  surf  zones  [e.g.,  Guza 
and  Thornton  (1980)].  The  oscillatory  components  of  the  horizontal  velocity  and 
free  surface  elevation  are  then  related  as  follows: 

{u-u)  ~  ^{rj -g)  (5.1) 

Here  g  =  gravitational  acceleration  and  h  =  mean  water  depth,  given  hy  h  =  {d  +  rj) 
with  d  =  still  water  depth.  If  equation  (5.1)  holds,  the  following  relationships  can 
be  derived  as  well: 


~  s 


(5.2) 


in  which  cr,  s  and  K  are  the  standard  deviation,  skewness,  and  kurtosis  of  the 
free  surface  elevation,  g.  Furthermore,  the  probability  density  functions  of  it*  = 
[u  —  u)  /au  and  »;*  =  {g  —  g)  fa  become  approximately  the  same  because  it*  ~  77* 
under  the  assumption  of  (5.1),  with  cr„  ~  [\J'gJh)a. 

To  estimate  the  undertow  velocity,  u,  the  horizontal  velocity,  u,  at  the  ap¬ 
proximate  mid-depth  is  assumed  to  be  represented  by  the  depth-averaged  velocity, 
U .  For  an  impermeable  beach,  the  time-averaged  continuity  equation,  hU  =  0, 
yields: 

U  =  -{hy'  (,-^)(t/-£7)  (5.3) 


in  which  h  =  {d  +  g)  (Kobayashi  et  al.  1989).  Substitution  of  u  ~  C/  and  equation 
(5.1)  into  (5.3)  yields: 


U  — 


(5.4) 


This  study  will  compare  predictions  of  equations  (5.2)  and  (5.4)  to  velocity  statistics 
obtained  from  tests  1  and  2. 


5.2  Cross-Shore  Variations  in  Statistical  Values  of  Horizontal  Velocity 

Figures  5.1  and  5.2  show  the  cross-shore  variations  of  the  measured  values 
of  u,  CTu,  Su  and  Ku  in  the  still  water  depth  d  =  15-58  cm  for  test  1  and  d  =  17- 
57  cm  for  test  2.  The  three  data  points  at  cross-shore  locations  x  =  7.35  m  and  x  = 
9.85  m  for  the  repeated  velocity  measurements  made  in  vertical  velocity  profiles  are 
represented  by  the  middle  point  at  each  location  to  be  consistent  with  the  other  data 
points  in  these  figures.  The  depth  range  was  limited  due  to  difficulties  measuring 
undertow  velocity,  u,  in  the  shallow  water.  The  computed  values  of  u,  (Tu,  Su,  and 
Ku  at  each  cross-shore  location  are  also  shown  in  Figures  5.1  and  5.2.  These  were 
obtained  from  equations  (5.2)  and  (5.4),  making  use  of  the  corresponding  measured 
values  of  /i  =  (d  -f  rj),  a,  s,  and  K. 

The  measured  and  computed  undertow,  u,  agree  well  in  deeper  water  in  both 
tests.  However,  there  is  some  underprediction  of  the  magnitude  of  u  in  the  surf 
zone.  This  may  be  explained  by  the  fact  that  equation  (5.4),  which  is  based  on 
linear  long  wave  theory,  does  not  account  for  the  additional  water  volume  flux  due 
to  a  roller.  This  roller  effect  increases  the  magnitude  of  u,  as  shown  by  Svendsen 
(1984)  for  regular  waves.  This  underprediction  was  also  visible  in  Tests  2  and  3  of 
Herrman  et  al.  (1997).  The  measured  undertow  velocity  fields  are  shown  over  the 
corresponding  profiles  in  Figure  5.3.  Undertow  is  strong  over  most  of  the  terrace 
region  in  both  tests,  leading  one  to  expect  offshore  transport.  This  will  be  examined 
further  in  Chapter  6. 

There  is  good  agreement  (within  about  10%)  between  the  measured  and 
computed  standard  deviation,  cr„,  for  test  1,  as  shown  in  Figure  5.1.  Figure  5.2 
shows  only  fair  agreement  (usually  within  30%)  in  test  2.  These  results  compare 
well  to  the  level  of  agreement  predicted  by  Guza  and  Thornton  (1980)  from  field 
data  intercomparisons.  It  should  be  noted,  however,  that  Guza  and  Thornton  (1980) 
focused  their  analysis  on  wind  wave  frequencies,  while  in  this  study  the  measured. 
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Figure  5.4:  Comparison  of  Linear  Theory  and  Finite  Depth  Theory  Predictions  of 
Standard  Deviation,  cr„,  for  Test  2 


unfiltered  time  series  are  used  because  of  difficulty  separating  wind  wave  and  low- 
frequency  components  for  practical  applications.  The  merely  fair  agreement  shown 
in  test  2  of  this  study  might  be  explained  by  examination  of  the  values  of  Lpjh, 
where  Lp  =  linear  wavelength  based  on  Tp  and  Tp  =  2.8  and  1.6  s  for  tests  1  and 
2,  respectively,  as  listed  in  Table  3.5.  For  depths  ranging  from  15  to  58  cm,  Lpjh 
ranged  from  11  to  22  in  test  1  but  only  from  6  to  12  in  test  2.  This  suggests  that 
the  linear  long-wave  assumption  may  not  be  appropriate  for  test  2.  Indeed,  when 
the  standard  deviation,  (t„,  is  calculated  for  test  2  using  linear  finite-depth  theory 
based  on  the  peak  period  Tp,  there  is  better  agreement  between  theory  and  data. 
This  is  illustrated  in  Figure  5.4. 

The  values  of  |u|  and  (7„  initially  increase  landward  in  Figures  5.1  and  5.2,  con¬ 
sistent  with  values  based  on  the  depth-averaged  velocity,  U,  computed  by  Kobayashi 
and  Karjadi  (1996)  using  the  finite-amplitude  shallow-water  equations.  However, 
unlike  the  1:16  slope  data  of  Herrman  et  al.  (1997),  both  |tt|  and  (Tu  decrease  land¬ 
ward  after  a  peak  near  the  location  of  the  terrace  edge.  The  intense  breaking  at 
the  bar  crest  in  test  1  appears  to  have  caused  the  sharp  decrease  of  |u|  and  cr^  in 
Figure  5.1.  In  summary.  Figures  5.1  and  5.2  indicate  that  the  local  application 


of  linear  long-wave  theory  is  a  reasonable  first  approximation  for  u  and  in  the 
shoaling  and  surf  zones. 

On  the  other  hand,  these  figures  generally  indicate  that  the  skewness  and 
kurtosis  (the  third  and  fourth  moments)  are  not  well  predicted  by  linear  long-wave 
theory.  The  skewness,  Su-,  of  the  horizontal  velocity,  u,  at  the  approximate  mid¬ 
depth  is  always  less  than  the  skewness,  s,  of  the  free  surface  elevation,  rj.  Also, 
drops  steeply  at  the  point  of  intense  breaking  in  test  1.  Milder  declines  in  due  to 
wave  breaking  were  measured  by  Herrman  et  al.  (1997)  in  their  tests  2  and  3.  The 
comparison  between  s  and  Su  suggests  that  nonlinearity  decreases  downward  from 
the  free  surface.  The  difference  (s  —  s„)  tends  to  fluctuate,  with  an  average  value  on 
the  order  of  0.5.  The  measured  values  of  Su  differ  noticeably  in  the  two  tests.  While 
Su  in  test  1  seems  to  have  been  drastically  affected  by  breaking  at  the  terrace  edge, 
the  shoreward  increase  of  s„  for  test  2  is  only  slightly  perturbed  by  breaking  over  the 
bar.  The  skewness,  s^,  in  the  shoaling  region  for  test  2  is  small  and  slightly  negative, 
which  implies  that  the  cross-shore  velocities  are  skewed  somewhat  seaward. 

The  kurtosis,  Ku,  of  the  mid-depth  horizontal  velocity,  u,  is  generally  over¬ 
predicted  by  the  kurtosis,  K,  of  the  free  surface  elevation,  rj,  in  both  tests  1  and 
2.  For  test  1,  the  largest  differences  between  the  data  and  the  linear  theory  occur 
in  the  shoaling  zone  in  deeper  water.  In  the  surf  zone,  kurtosis  values  show  greater 
variability  for  test  1,  again  reflecting  the  intense  breaking  that  took  place  in  that 
test.  For  test  2,  Ku  is  less  than  K  at  all  locations,  with  only  a  small  shift  in  values 
to  mark  the  surf  zone.  Both  tests  here  indicate  closer  agreement  than  the  kurtosis 
comparison  by  Herrman  et  al.  (1997). 

The  values  plotted  in  Figures  5.1  and  5.2  are  also  listed  in  Tables  5.1  and  5.2, 
where  Xu  and  Zu  are  the  horizontal  and  vertical  coordinates  of  the  ADV  probe.  The 
coordinate  system  used  herein  has  been  defined  in  Figure  3.2. 


Table  5.1:  Velocity  Data  for  Test  1 
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1-20  -18.3  10.35  -0.14  -7.74  -5.41  21.94  26.32  0.29  0.86  3.26  3.35 

1-21  -15.5  10.85  -0.14  -6.80  -6.50  23.84  27.82  0.19  0.89  3.02  3.64 


Table  5.2:  Velocity  Data  for  Test  2 


^  o 

CS  CO 
CO  CO 
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-3.27  -  16.55  -  0.34  —  2.45 

-4.28  -3.90  17.13  22.52  0.41  0.96  2.50 


5.3  Vertical  Variations  of  Velocity  Statistics  at  Two  Cross-Shore 
Locations 

Vertical  profiles  of  the  cross-shore  velocity  were  taken  at  a:  =  7.35  m  and 
9.85  m  to  examine  the  variation  of  «,  (T„,  and  over  depth.  At  a:  =  7.35  m 
velocity  time  series  were  recorded  at  elevations  Zu  =  -12,  -15,  and  -18  cm  below 
SWL  for  both  tests.  At  a:  =  9.85  m  data  was  taken  at  =  -8,  -14,  and  -19  cm 
below  SWL  for  test  1  and  =  -8,  -14,  and  -18  cm  below  SWL  for  test  2.  These 
results  have  already  been  listed  in  Tables  5.1  and  5.2.  In  this  section  they  will  be 
presented  graphically.  Figures  5.5  and  5.6  show  profiles  of  u,  (T„,  s„,  and  Ku  at  the 
two  locations  for  tests  1  and  2,  respectively.  The  locations  of  these  velocity  profiles 
with  respect  to  the  equilibrium  beach  profiles  may  be  viewed  in  Figure  5.3. 

In  general,  the  results  indicate  that  the  mean  velocity  and  standard  deviation 
remain  relatively  constant  over  the  depth  ranges  measured.  Measured  undertow, 
u,  is  within  20%  of  its  middepth  value  for  all  profiles.  The  vertical  variation  of 
undertow  u  for  irregular  waves  appears  to  be  less  than  that  for  regular  waves  (Cox 
and  Kobayashi  1997).  Standard  deviation,  <Tu,  remains  within  5%  of  the  value 
calculated  at  middepth  for  both  tests. 

Figures  5.5  and  5.6  also  show  relatively  uniform  values  for  the  higher  mo¬ 
ments,  Su  and  Ku,  over  the  depth.  There  is  greater  variability  in  these  results  at 
X  =  9.85  m  for  test  1,  however,  which  might  be  explained  as  a  consequence  of  the 
intense  wave  breaking.  Variability  may  also  be  due  to  instrument  error  in  measuring 
the  velocity,  especially  at  the  shallowest  depth,  Zu  =  -8  cm.  Except  for  this  one  po¬ 
sition,  skewness  Su  varied  less  than  4  %  from  its  middepth  value,  and  kurtosis,  Ku, 
remained  within  3%  of  its  middepth  value  in  all  other  tests.  The  vertical  profiles 
of  cross-shore  velocity  and  its  moments  are  important  in  predicting  the  cross-shore 
sediment  transport  using  the  measured  velocity  at  a  single  elevation,  as  will  be 
discussed  in  Chapter  6. 


Depth  (cm)  Depth  (cm) 
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Figure  5.5:  Velocity  Profile  Statistics  for  Test  1 
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Figure  5.6:  Velocity  Profile  Statistics  for  Test  2 


5.4  Measured  and  Theoretical  Probability  Density  Functions 

In  this  section  the  measured  probability  density  function,  /(u»),  is  compared 
with  the  exponential  gamma  distribution  based  on  the  measured  skewness,  The 
agreement  is  good  in  test  1  except  for  slight  deviations  at  the  peaks  of  /(u*);  these 
deviations  are  somewhat  larger  in  the  results  for  test  2.  To  show  the  degree  of  the 
overall  agreement  for  tests  1  and  2  together,  the  combined  measured  distributions 
are  separated  into  two  groups  on  the  basis  of  the  measured  skewness,  Figure  5.7 
shows  the  measured  distributions  with  =  (-0.14)  -  0.29  in  comparison  with  the 
exponential  gamma  distributions  with  =  0  and  =  0.29,  where  >  0  must  be 
assumed  to  be  consistent  with  the  assumption  of  s  >  0  made  in  (2.1).  Data  falling 
outside  the  theoretical  curves  most  often  corresponded  to  values  of  <  0.  In  the 
lower  plot,  the  measured  distributions  with  s„  =  0.31  -  0.62  are  compared  with  the 
exponential  gamma  distributions  with  =  0.31  and  0.62.  These  exponential  gamma 
distributions  deviate  somewhat  more  within  the  specified  range  of  Figure  5.7 
indicates  that  most  of  the  data  points  fall  between  these  two  theoretical  curves, 
apart  from  the  scatter  near  the  peak  of  /(u*).  The  agreement  is  similar  to  those 
shown  in  Figures  4.4  and  4.5  for  the  free  surface  elevations.  Lack  of  sufficient  data 
precludes  determining  whether  the  exponential  gamma  distribution  is  applicable  to 
horizontal  velocities  near  the  still  water  shoreline.  For  completeness,  the  full  set  of 
individual  comparisons  is  shown  in  Section  5.6  in  Figures  5.11  -  5.22. 

5.5  Synchronization  of  Velocity  and  Free  Surface  Time  Series 

As  depicted  in  Figure  3.2,  the  acoustic  Doppler  velocimeter  (ADV)  and  wave 
gauge  5  were  positioned  at  the  same  cross-shore  location  for  each  run  of  each  test 
in  order  to  measure  simultaneous  velocity  and  free  surface  time  series  at  each  point. 
Such  a  synchronization  would  enable  the  separation  of  incident  and  reflected  waves 
in  shallow  water  without  using  multiple  wave  gauges.  Unfortunately,  attempts  at 
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Figure  5.7:  Measured  and  Computed  Probability  Distributions  of  Horizontal  Ve¬ 
locities  for  Tests  1  and  2  in  Two  Separate  Skewness  Ranges 


synchronization  were  ultimately  unsuccessful.  Proper  synchronization  was  ham¬ 
pered  by  the  need  for  two  separate  computers  to  run  the  ADV  and  the  wave  gauges. 
Different  clock  and  program  speeds  on  the  two  computers  introduced  slight  phase 
shifts  into  the  data  so  that  after  400  s  of  testing,  velocity  and  free  surface  time 
series  that  were  initially  synchronized  had  shifted  by  approximately  0.33  s.  The 
development  of  such  a  phase  shift  is  illustrated  in  Figure  5.8. 

The  figure  shows  the  400  s  time  series  of  the  normalized  free  surface  elevation, 
7/*  =  (7/  —  »f)/cr,  and  the  normalized  cross-shore  velocity  u*  =  (u  —  u)/(T„,  for  the 
first  run  of  test  1  in  three  segments,  where  rf,  cr,  u,  and  are  the  measured  values 
for  this  run.  The  first  25  s  of  both  time  series  are  presented  in  expanded  form  in  the 
top  panel  to  illustrate  the  close  synchronization  of  the  initial  measurements.  In  the 
second  panel,  the  next  325  s  of  the  two  time  series  are  displayed  in  a  compacted  form. 
Careful  examination  of  this  panel  reveals  increasing  evidence  of  the  phase  shift  as 
one  moves  from  left  to  right.  The  third  and  final  panel  displays  the  remaining  25  s  of 
the  time  series  in  expanded  form  to  show  clearly  the  phase  shift  that  has  developed. 
Figure  5.9  shows  a  plot  of  the  individual  measurement  points  of  us.those  of  7/, 
for  the  same  run  as  above.  The  general  trend  of  the  data  points  follows  the  line 
«*  =  77*,  corresponding  to  equation  (5.1)  based  on  local  linear  progressive  long- wave 
theory. 

It  was  initially  suspected  that  the  phase  shift  may  have  been  due  primarily  to 
the  effects  of  reflected  waves.  However,  more  careful  testing  revealed  the  computer 
speeds  to  be  almost  entirely  responsible  for  the  time  series  difference.  A  more 
accurate  measure  of  the  phase  shift  was  obtained  by  running  the  ADV  and  wave 
gauge  data  programs  simultaneously  for  3  hours.  After  this  time  the  ADV  computer 
time  reading  was  lagging  behind  that  of  the  wave  gauge  computer  by  9  seconds.  The 
phase  shift  rate  implied  by  this  time  difference,  0.33  s  per  400  s,  was  used  to  adjust 
time  series  of  77*,  as  illustrated  in  Figure  5.8;  given  this  adjustment,  very  close 
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sychronicity  of  the  time  series  was  obtained.  The  corrected  version  of  the  above 
time  series  of  ?/*  is  replotted  with  the  original  time  series  of  tt*  in  Figure  5.10.  The 
final  25  s  of  the  time  series  still  do  not  show  a  perfect  match  between  u*  and  77*. 
These  additional  irregularities  are  most  likely  due  to  the  reflected  waves  which  had 
reached  the  gauge  position  by  that  time  but  were  not  present  for  the  initial  25  s  of 
measurements. 

It  would  be  possible  to  correct  the  actual  time  series  data  using  the  phase 
shift  determined  above.  This  would  require  the  use  of  a  point-by-point  interpolation 
program  to  shift  each  of  the  8000  points  of  the  7;*  or  tt*  time  series  by  the  appro¬ 
priate  amount,  based  on  its  location  in  the  series.  However,  the  validity  of  such  an 
adjustment  would  be  questionable.  It  might  not  be  possible  to  verify  the  accuracy 
of  the  estimated  0.33  s  shift,  and  any  errors  in  this  measurement  would  propagate 
through  the  analysis.  For  these  reasons  it  was  decided  not  to  continue  with  efforts 
to  synchronize  the  free  surface  and  velocity  time  series. 


395 


400 


390 

rime  (sec) 

me  Series  of  Normalized  Free  Surface,  7;», 
t*,  Run  1  of  Test  1 


3 


5.6  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity 

In  Figures  5.11  -  5.22  the  probability  distributions  of  the  horizontal  velocity 
are  compared  with  the  distributions  computed  from  the  exponential  gamma  distri¬ 
bution  at  each  of  the  21  velocity  measurement  locations  for  tests  1  and  2.  Note  that 
three  separate  distributions  corresponding  to  different  elevations  are  shown  for  the 
locations  x  =  7.35  m  and  x  -  9.85  m,  where  velocity  profiles  were  measured.  The 
probe  coordinates,  and  z^,  used  in  Tables  5.1  and  5.2  are  indicated  as  xu  and  zu 
in  each  diagram. 


Probability  Density  Function  f 


Figure  5.11:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Positions  x  =  2.85,  3.35,  3.85,  and  4.35  m  for  Test  1 
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Probability  Density  Function  f 


Figure  5.12:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Positions  x  =  4.85,  5.35,  5.85,  and  6.35  m  for  Test  1 
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Figure  5.13:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Positions  x  =  6.85  and  7.35  m  for  Test  1  (Measurements 
Repeated  Due  to  Three-Point  Velocity  Profile  at  7.35  m) 


Probability  Density  Function  f 


Figure  5.14:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Positions  x  =  7.85,  8.35,  8.85,  and  9.35  m  for  Test  1 
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Test  1 


Figure  5.15:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Positions  x  =  9.85  and  10.35  m  for  Test  1  (Measurements 
Repeated  Due  to  Three-Point  Velocity  Profile  at  9.85  m) 


Test  1 


Figure  5.16;  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Position  x  =  10.85  m  for  Test  1 


Test  2 


Figure  5.17:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Positions  x  =  2.85,  3.35,  3.85,  and  4.35  m  for  Test  2 


Test  2 


Figure  5.18:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Positions  x  =  4.85,  5.35,  5.85,  and  6.35  m  for  Test  2 


Test  2 


Figure  5.19:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Positions  x  =  6.85  and  7.35  m  for  Test  2  (Measurements 
Repeated  Due  to  Three-Point  Velocity  Profile  at  7.35  m) 
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5.20:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Positions  x  =  7.85,  8.35,  8.85,  and  9.35  m  for  Test  2 
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Figure  5.21:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Positions  x  =  9.85  and  10.35m  for  Test  2  (Measurements 
Repeated  Due  to  Three-Point  Velocity  Profile  at  9.85  m) 


Test  2 


Figure  5.22:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Position  x  =  10.85  m  for  Test  2 


Chapter  6 


CROSS-SHORE  SEDIMENT  TRANSPORT 

As  mentioned  in  Chapter  1,  perturbations  of  equilibrium  profiles  have  been 
observed  in  the  field  and  the  laboratory,  but  the  processes  involved  in  the  formation 
of  these  perturbed  beach  profiles  are  still  poorly  understood.  This  chapter  will  ex¬ 
amine  the  sediment  transport  models  based  on  the  sheet  flow  approach  [Trowbridge 
and  Young  (1989)]  and  the  energetics  approach  [Bowen  (1980);  Bailard  (1981)].  Af¬ 
ter  a  summary  of  theoretical  background,  the  necessary  parameters  for  each  model 
will  be  calculated  using  the  data  from  tests  1  and  2.  It  will  be  shown  that  the 
energetics-based  model  predicts  the  equilibrium  profile  better  than  the  sheet  flow 
model  for  both  tests.  Both  models  overpredict  the  rate  of  change  of  bottom  ele¬ 
vation,  dzbjdt,  in  the  region  of  the  bar  crest  for  test  1,  but  the  energetics  model 
predictions  are  within  1  cm/hr  at  all  locations  for  test  2.  The  shortcomings  of  these 
two  approaches  will  be  used  to  explore  the  reasons  why  existing  models  cannot 
predict  zero  net  sediment  transport  rate  on  these  equilibrium  profiles. 

6.1  Sheet  Flow  Model  for  Onshore  Bar  Movement 

The  sheet  flow  model,  developed  by  Trowbridge  and  Young  (1989),  is  presently 
the  only  model  available  that  attempts  to  explain  onshore  bar  movement  outside 
the  surf  zone.  For  the  profiles  measured  in  this  study,  there  may  be  some  regions 
where  the  sheet  flow  model  does  not  apply.  While  this  model  is  generally  limited  to 
smooth  beds  without  ripples,  the  profiles  obtained  in  these  experiments  did  include 


rippled  sand  beds  in  the  offshore  and  surf  zones  for  both  tests.  However,  both  pro¬ 
files  were  free  of  ripples  in  the  region  over  the  bar.  Hence,  it  is  likely  that  at  least 
quasi-sheet  flow  conditions  were  achieved  in  the  regions  of  intense  wave  breaking. 
Consequently,  the  following  comparison  will  indicate  whether  this  sheet  flow  model 
is  applicable  in  the  rippled  region  outside  the  surf  zone  as  well  as  in  the  sheet  flow 
and  rippled  regions  inside  the  surf  zone. 


6.1.1  Theory 

The  measured  cross-shore  velocity  is  represented  here  by  u{t),  and  the  time- 
averaged  cross-shore  velocity  or  undertow  is  represented  by  u.  The  instantaneous 
rate  of  onshore  sediment  transport,  q{t),  is  assumed  to  be  expressible  as 

(6.1) 

wd  pgd{s  —  1) 

where  w  =  fall  velocity,  d  =  diameter  of  uniform  sediment  which  is  taken  as  d  =  dso, 
p  =  density  of  fluid,  s  =  ps! p  —  specific  gravity  of  sediment  with  ps  =  sediment 
density,  n{t)  =  instantaneous  bottom  shear  stress,  g  =  gravitational  acceleration, 
and  K  =  &n  empirical  coefficient.  The  time-averaged  cross-shore  transport  rate, 
taken  to  be  positive  onshore,  may  thus  be  written  as 


q:=K- 


WTb 


(6.2) 
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Trowbridge  and  Young  (1989)  analyzed  the  wave  boundary  layer  and  derived  the 
following  expression  for  the  mean  bottom  shear  stress,  1%: 


—  ,  pfw 

Ti,  =  H - 


(6.3) 


where  f-u,  =  friction  coefficient  and  u  =  oscillatory  part  of  first-order  wave  velocity 
—  {u  —  u).  As  this  expression  only  includes  |u|,  it  does  not  explicitly  account  for 
possible  undertow,  which  was  significant  in  the  present  experiment.  Note  that  while 


Trowbridge  and  Young  defined  the  x-axis  as  positive  offshore,  the  present  analysis 
defines  x  as  positive  onshore.  Substituting  (6.3)  into  (6.2): 

'KU 


\  2  )  g{s-  l)^f^ 

Trowbridge  and  Young  calibrated  the  value  of  Kfu,  for  d  ~  0.2  mm  and  recom¬ 
mended  Kfu,  =  0.5.  In  the  present  experiment,  dso  =  0.18  mm  as  explained  in 
Section  3.1.1,  and  Kf^  =  0.5  will  be  adopted  as  well.  For  sheet  flow  conditions, 
equation  (6.4)  may  be  used  to  predict  the  time-averaged  cross-shore  transport  rate, 
q^,  based  on  velocity  data  as  indicated  by  the  subscript  u. 

Trowbridge  and  Young  used  linear  long  wave  theory  to  relate  u  to 


u  =  (u  -  u)  ~  y-(7/  -r/)  = 


(6.5) 


with  the  assumption  that  fj  is  Gaussian.  Since  velocity  data  is  available  here,  it  will 
be  assumed  instead  that  u  is  Gaussian  and  expressible  as 


m 


-(|uP)^-^ 


(6.6) 


This  assumption  is  better  than  that  of  Trowbridge  and  Young  because  the  skewness, 
of  the  velocity  is  generally  smaller  than  the  skewness,  s,  of  the  free  surface. 
Equation  (6.5)  then  implies 

(6.7) 


which  can  be  rewritten  as 


up  ~  ^|7)P 

^  9  ^2 

U 

Combining  (6.6)  and  (6.7),  one  obtains 


i^) 


(6.8) 


(6.9) 


This  expression  for  |uP  may  be  substituted  into  equation  (6.4)  to  obtain,  after  some 
algebra: 

/o 

(6.10) 
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Given  ®  =  cr^  and  a  =  i/^rms/v^,  as  in  equation  (2.8),  one  finally 

obtains  ^ 

_  K fm  f Hxias'^  /'n  1  1  ^ 

“  16v^(s-l)  \1~ )  ^  ’ 

Equation  (6.11)  may  be  used  to  predict  the  time- averaged  cross-shore  transport  rate, 

using  free  surface  data  as  indicated  by  the  subscript  q. 

An  expression  for  the  predicted  rate  of  change  of  the  sand  bottom  elevation 

may  be  derived  by  applying  the  conservation  equation  of  sediment  to  a  small  sand 

volume  of  depth  Azi,  length  Ax,  and  unit  width.  The  increase  in  sand  volume  over 

time  At  is  equal  to  the  net  influx  of  sand  during  that  time: 


(1  -  np)Azi,Ax  = 


At 


(6.12) 


In  the  limit  as  At  0,  (6.12)  may  be  rearranged  to  express  the  erosion  (negative)  or 
accretion  (positive)  rate  dzijdt  in  terms  of  the  gradient  of  the  cross-shore  transport 
rate  q\ 

^  ^  (6.13) 

dt  (1  —  rip)  dx 

in  which  rip  =  porosity.  Equation  (6.13)  may  be  combined  with  equation  (6.4)  or 
(6.11)  to  determine  the  erosion/accretion  rate  based  on  velocity  data  or  free  surface 
data,  respectively,  where  [dzbfdt)u  and  (dzb/dt)p  are  used  for  the  computed  values 
of  [dzbjdt)  corresponding  to  and  q^,  respectively. 


6.1.2  Analysis  of  Experimental  Results 

The  velocity  and  free  surface  data  from  tests  1  and  2  are  analyzed  to  obtain 
the  quantities  necessary  for  the  evaluation  of  equations  (6.4)  and  (6.11)  at  the  17 
cross-shore  locations  at  which  velocity  time  series  were  measured.  At  the  locations 
of  the  three  point  velocity  profiles  the  time  series  from  only  the  middle  position  are 
used;  the  vertical  variations  will  be  discussed  in  Section  6.2.2.  This  analysis  uses  the 
results  of  sand  tests  summarized  in  Section  3.1:  w  =  1.9  cm/s,  d  =  d^o  =  0.18  mm. 


Figure  6.1:  Ratio  (|uP)/(|«P)^’^  for  Tests  1  and  2 


s  =  2.66,  and  Up  =  0.4.  A  cubic  spline  interpolation  is  performed  on  each  set  of  q 
values  to  obtain  the  corresponding  derivatives,  dqldx.  Equation  (6.13)  is  then  used 
to  calculate  dz\,ldt  at  each  location.  The  results  of  these  calculations  are  displayed 
graphically  in  Figures  6.1  -  6.3  and  summarized  in  Tables  6.1  and  6.2  at  the  end  of 
this  section. 

Figure  6.1  provides  a  check  on  equation  (6.6),  which  was  based  on  the  as¬ 
sumption  that  u  is  Gaussian.  This  equation  may  be  rewritten  as 


u 


f8 


=  G-~1.60 
(|u|2)l-5  V  TT 


(6.14) 


Figure  6.1  shows  that  the  ratio  is  indeed  quite  close  to  this  Gaussian  prediction  for 
both  tests  at  nearly  all  locations,  with  some  minor  deviation  in  test  1  in  the  region 
of  broken  waves. 

As  mentioned  above,  the  preconditions  of  the  sheet  flow  model  may  be  sat¬ 
isfied  over  only  part  of  these  equilibrium  profiles  because  of  the  ripples  that  were 
present.  The  sheet  flow  conditions  occurred  from  about  a;  =  5mtox  =  9m. 

Figures  6.2  and  6.3,  along  with  Tables  6.1  and  6.2,  predict  large  onshore 
transport  rates  and  rapid  profile  change  in  the  surf  zone,  especially  for  test  1.  For 
perfectly  equilibrium  profiles,  dzhldt  =  0,  but  the  measured  equilibrium  profiles  for 
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Figure  6.3:  Net  Transport  Predictions  and  Bottom  Elevation  Change  Rates 

(dzbldt)u,  {dzbfdt)j„  Over  Test  2  Profile  (Sheet  Flow  Model) 
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test  1  and  2  discussed  in  Section  3.4  had  uncertainties  on  the  order  of  1  cm/hr  or 
less.  For  test  2,  dzijdt  based  on  u  is  somewhat  smaller  than  dzhjdt  based  on  rj, 
though  still  greater  than  1.0  cm/hr  over  the  bar.  For  both  tests  the  bar  is  predicted 
to  move  further  shoreward.  Rapid  change  of  this  sort  was  of  course  not  observed  in 
any  part  of  either  profile.  However,  it  should  be  stated  that  this  sheet  flow  model 
actually  predicts  the  values  of  dzhjdt  well  below  the  data  uncertainty  of  1  cm/hr 
outside  the  surf  zone,  even  though  ripples  were  present  outside  the  surf  zone. 

To  apply  the  model  of  Trowbridge  and  Young  (1989)  in  the  surf  zone,  the 
effects  of  wave  breaking  and  undertow  will  need  to  be  included.  These  played  an 
important  part  in  surf  zone  dynamics  in  the  experiment  and  were  illustrated  in 
Figure  5.3.  The  inclusion  of  such  effects  would  act  to  oppose  the  purely  wave- 
induced  onshore  transport  of  this  model  and  would  result  in  the  prediction  of  a 
more  equilibrium  profile. 


6.2  Energetics-Based  Model  for  Offshore  Bar  Movement 

Unlike  the  sheet  flow  model,  the  energetics-based  model  developed  by  Bowen 
(1980)  and  Bailard  (1981)  attempts  to  account  for  both  onshore  transport  due  to 
wave  skewness  and  offshore  transport  due  to  undertow,  as  well  as  slope  effects  due  to 
gravity.  This  model  separates  sediment  transport  into  bed  load  and  suspended  load 
components,  including  separate  terms  for  each  of  the  above  transport  mechanisms. 
Most  recently,  Thornton  et  al.  (1996)  used  the  energetics-based  model  to  predict  bar 
movement  in  the  field  at  Duck,  North  Carolina.  Sediment  size  ranged  from  0.12- 
0.20  mm,  and  significant  wave  height  reached  2. 0-2. 5  m.  Their  results  indicated  that 
the  model  could  predict  correctly  the  bar  movement  during  storms  to  first  order,  but 
that  it  underpredicted  the  development  of  a  trough  and  performed  poorly  in  mild 
wave  conditions.  They  found  that  the  effects  of  bottom  slope  on  sediment  load  were 
negligible,  while  the  suspended  load  was  calculated  to  be  an  order  of  magnitude 
greater  than  the  bed  load. 

In  the  present  experiment,  the  fine  sand  (dso  =  0.18  mm)  was  approximately 
the  same  size  as  that  present  at  Duck,  while  the  waves  {Hs  ^  0.2  m)  were  an  order 
of  magnitude  smaller,  though  still  relatively  large  for  laboratory  tests.  The  terraced 
and  barred  beaches  for  tests  1  and  2  were  similar  to  the  beach  profiles  measured 
in  Duck,  except  that  the  offshore  slope  in  these  tests  was  much  steeper.  Because 
the  experiments  were  conducted  in  a  wave  flume  there  was  no  longshore  current  to 
contribute  to  sediment  transport.  The  following  analysis  will  refrain  from  separating 
the  oscillatory  velocity  into  high  and  low  frequency  components  and  will  examine 
only  the  unfiltered  velocity  time  series. 

6.2.1  Theory 

In  the  energetics  model,  the  time-averaged  cross-shore  sediment  transport 
rate  per  unit  width  can  generally  be  expressed  by  equation  (2)  in  Thornton  et  al. 


(1996).  For  the  present  analysis  with  zero  longshore  velocity,  however,  the  net 
onshore  sediment  transport  rate  q  is  given  by: 

q  =  (|M(<)|*S(i))  +  ATi  (|u(()Pu)  -  ATij  (|M(i)P) 

+  K,  (Rijpiw)  +  K.  (Ri)p“)  -  K.,  (|ti(()l')  (6.15) 

—  Qbw  ^bu  Qsw  d”  Qsu  d”  Qsg 

where  u{t)  =  cross-shore  horizontal  velocity,  which  is  positive  shoreward  in  this 
analysis.  In  equation  (6.15),  the  first  three  terms  represent  bed  load  (subscript 
b)  produced  by  wave  asymmetry  (subscript  w),  undertow  (subscript  u),  and  the 
effects  of  gravity  (subscript  g)  on  the  bottom  slope,  respectively.  The  final  three 
terms  represent  suspended  load  (subscript  s)  produced  by  the  same  three  respective 
effects.  The  coefficients  in  equation  (6.15)  are  expressed  as 

-  -J—r  -  K 

^  (s  —  l)p  ^ *’tan((^) 

Ks  =  ;  Ks,=Ks-tM^) 

{s-l)g  w  w 

where  s  =  ps/p  =  specific  gravity  of  sand,  Cf  =  bed  drag  coefficient,  (f)  =  internal 
friction  angle  of  sand,  Cb  =  bed  load  efficiency  factor,  £«  =  suspended  load  efficiency 
factor,  tan(;S)  =  local  bed  slope  =  dzbjdx,  and  w  =  fall  velocity.  In  the  following 
analysis,  s  =  2.66  and  w  =  1.9  cm/s,  and  other  parameters  are  given  the  same 
values  as  those  used  in  Thornton  et  al.  (1996):  Cf  —  0.003,  tan((^)  =  0.63,  £&  = 
0.135,  and  =  0.015.  The  local  slope,  tan()S),  is  computed  using  the  equilibrium 
bottom  profile,  Zb{x),  for  each  test. 

After  q  is  obtained  using  equation  (6.15),  the  predicted  erosion  or  accretion 
rate,  dzhjdt,  for  the  profile  at  given  locations  is  found  by  use  of  equation  (6.13). 
This  corresponds  to  equation  (1)  in  Thornton  et  al.,  except  that  here  the  value  of 
^  =  (1  —  Up)  is  taken  as  0.6  rather  than  0.7  because  the  measured  porosity  Up  =  0.4 
in  these  experiments. 


Table  6.3;  Mean  Velocity  Statistics,  Test  1 
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6.2.2  Analysis  of  Experimental  Results 

Each  of  the  time-averaged  velocity  expressions  in  equation  (6.15)  was  eval¬ 
uated  for  the  21  velocity  measurements  where  the  middle  point  at  the  locations  of 
the  three-point  velocity  profiles  is  used  for  plotting  the  cross-shore  variations.  These 
values  are  tabulated  for  tests  1  and  2  in  Tables  6.3  and  6.4.  The  values  of  the  six 
sediment  load  components  are  plotted  over  the  test  1  and  test  2  equilibrium  profiles 
in  Figures  6.4  and  6.5  and  listed  in  Tables  6.5  and  6.6.  In  both  tests  the  largest 


Table  6.4:  Mean  Velocity  Statistics,  Test  2 
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values  predicted  for  all  types  of  sediment  transport  occur  around  the  location  of 
the  bar  or  terrace  edge,  in  the  vicinity  of  x  ~  7-8  m.  The  more  intense  breaking 
of  test  1  is  clearly  visible  from  a  comparison  of  sediment  loads  in  this  region.  The 
effects  of  bottom  slope  on  sediment  transport  {q^g,  %g)  are  most  clear  seaward  of 
the  bar,  where  the  bottom  slope  is  relatively  steep  and  positive.  Both  bed  and  sus¬ 
pended  sediment  fluxes  due  to  slope  effects  are  consistently  negative  (z.e.,  offshore 
transport).  Shoreward  of  the  bar,  where  the  bottom  slope  is  negative  for  a  short 
distance,  sediment  fluxes  due  to  slope  effects  are  noticeably  positive,  especially  in 
test  2.  One  additional  effect  illustrated  by  Figures  6.4  and  6.5  is  also  noteworthy. 
hi  X  —  7.85  m  in  test  1,  a  large  volume  of  sediment  is  transported  seaward  by  both 
wave  and  undertow  effects;  however,  just  0.5  m  shoreward  at  x  =  8.35  m,  the  wave 
suspended  load  suddenly  reverses  and  becomes  strongly  positive  onshore. 

Tables  6.5  and  6.6  list  three  sets  of  values  calculated  at  the  locations  x  = 
7.35  m  and  x  =  9.85  m,  where  three-point  velocity  profiles  were  taken.  The  predicted 
total  bed  load,  q^  =  -f-  q^^  +  qbg)^  and  the  predicted  total  suspended  load, 

-|-  q^^  -f-  q^g),  are  also  tabulated.  In  general,  the  predicted  suspended  load 
values  are  larger  than  the  corresponding  bed  load  terms,  as  would  be  expected  from 
the  relatively  large  waves  and  fine  sand  used  in  these  experiments.  However,  unlike 
Thornton  et  al.  (1996),  suspended  load,  q^,  is  rarely  an  entire  order  of  magnitude 
larger  than  bed  load,  in  these  laboratory  experiments. 

The  largest  suspended  sediment  fluxes  computed  in  these  experiments  are 
about  0.4  m^/day,  while  those  obtained  by  Thornton  et  al.  approached  50  m^/day. 
The  largest  bed  load  fluxes  in  these  experiments  were  approximately  0.1  m^/day, 
while  the  field  experiments  showed  values  no  larger  than  4  m^/ day.  Large  amounts 
of  sediment  were  predicted  to  have  been  suspended  by  a  strong  longshore  current  in 
the  field  experiments.  Also,  because  of  the  relatively  steeper  offshore  slope  in  the 
laboratory,  the  effects  of  bottom  slope  were  probably  proportionately  greater  in  the 
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Table  6.5:  Partial  Sediment  Loads  With  Energetics  Model,  Test  1 
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Table  6.6:  Partial  Sediment  Loads  With  Energetics  Model,  Test  2 
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lab  than  in  the  field.  In  tests  1  and  2  the  bed  load  due  to  slope  effects,  is  larger 
than  the  suspended  load,  q^g,  for  most  locations.  Moreover,  the  total  bed  load  qi,  is 
actually  comparable  to  the  total  suspended  load,  q^,  at  most  test  2  locations.  This 
may  be  reasonable  under  the  milder  wave  conditions  of  the  second  test. 

As  velocity  time  series  were  measured  at  three  elevations  at  two  cross- shore 
locations  in  these  experiments,  it  is  also  possible  to  examine  whether  the  cross¬ 
shore  sediment  transport  rate  q  in  equation  (6.15)  is  sensitive  to  the  elevation  of 
the  cross-shore  velocity,  u(f),  used  for  its  prediction.  This  variation  is  illustrated  in 
Figures  6.6  -  6.9.  While  predicted  sediment  loads  due  to  undertow  and  bottom  slope 
remain  relatively  constant  over  the  depth,  predictions  of  wave-induced  loads  q^^  and 
qf^^  vary  somewhat,  especially  at  x  =  9.85  m  for  test  1.  It  may  thus  be  possible  to 
predict  sediment  load  due  to  undertow  and  bottom  effects  by  measuring  velocities 
at  mid-depth  instead  of  immediately  outside  the  bottom  boundary  layer  as  specified 
by  the  theory  of  Bailard  (1981).  However,  for  predictions  of  wave-induced  load  q^^ 
and  qi,^,  the  use  of  mid-depth  velocity  measurements  may  cause  some  errors  when 
compared  to  velocity  measurements  made  near  the  bottom. 

The  total  sediment  loads  predicted  by  the  energetics  model  at  each  cross¬ 
shore  location  are  plotted  over  the  equilibrium  profiles  in  Figures  6.10  and  6.11. 
The  values  in  these  figures  are  summarized  in  Tables  6.7  and  6.8,  along  with  the 
ratio  of  total  suspended  load  to  total  bed  load,  and  the  predicted  erosion 

or  accretion  rate,  dzbfdt.  As  indicated  earlier,  the  ratio  q^jqb  rarely  reaches  the 
order  of  10  in  test  1  and  remains  fairly  close  to  the  order  of  unity  for  all  locations 
in  test  2.  This  may  be  reasonable  compared  to  the  field  data  of  Thornton  et  al. 
(1996)  because  wave  heights  in  the  laboratory  were  much  smaller.  In  test  1,  the 
net  sediment  transport,  q,  is  predicted  to  be  weakly  onshore  at  locations  seaward 
of  the  breaker  zone  but  strongly  offshore  at  almost  all  locations  on  the  terrace. 
The  magnitude  of  the  parameter  dzbjdt  is  not  greater  than  1  cm/hr,  except  in  the 
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Figure  6.6:  Variation  of  Suspended  Load  and  Bed  Load  Quantities  Over  Depth  at 
X  =  7.35  m,  Test  1  (Energetics  Model) 


region  of  intense  breaking.  In  test  2,  all  locations  are  predicted  to  have  a  weak 
offshore  transport,  with  dzbfdt  consistently  small.  These  values  are  smaller  than 
those  in  test  1  because  waves  broke  less  intensely,  leading  to  slower  profile  changes. 
Energetics  theory  thus  may  be  said  to  be  capable  of  predicting  equilibrium  in  test 
2  within  measurement  errors  of  1  cm/hr. 

For  both  tests  the  energetics  model  predicts  the  growth  of  a  bar  near  x 
=  7  m  and  a  trough  near  a;  =  8  m  at  rates  greater  than  those  observed  in  the 
laboratory.  It  should  be  noted  that  Bailard  (1981)  limits  the  applicability  of  his 
model  to  a  plane,  sloping  beach  and  does  not  mention  the  presence  of  ripples  or 
bars.  The  dominant  mechanism  for  the  development  and  movement  of  the  bar  in 
measurements  conducted  by  Thornton  et  al.  (1996)  was  strong  longshore  current 
and  undertow,  yet  in  the  experiments  it  has  been  shown  that  barred  and  terraced 
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Figure  6.7:  Variation  of  Suspended  Load  and  Bed  Load  Quantities  Over  Depth  at 
X  =  9.85  m,  Test  1  (Energetics  Model) 

profiles  may  develop  and  stabilize  without  longshore  current  but  with  relatively 
large  undertow.  While  the  energetics  model  underpredicted  trough  development  for 
Thornton  et  al.,  it  overpredicts  both  the  trough  and  bar  development  for  the  tests 
in  this  study. 

The  energetics  model  is  more  accurate  than  the  sheet  flow  model  in  both  tests 
1  and  2,  largely  because  of  the  inclusion  of  a  greater  number  of  sediment  transport 
mechanisms  in  the  former  model.  However,  both  models  assume  the  instantaneous 
response  of  both  bed  load  and  suspended  load  particles  to  the  horizontal  velocity 
immediately  outside  the  bottom  boundary  layer.  The  water  was  observed  to  be 
consistently  very  cloudy  during  both  tests,  and  the  assumption  of  the  instantaneous 
response  of  suspended  load  appears  to  be  questionable. 
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Test  2,  Xu=7.35m 


-0.05 

1 

1 

1 

1 

1 

1 

1 

1 

i 

1 

1 

E* 

1 

1 

•  ' 

1 

1 

1 

•  ' 

1 

1 

L 

-0.15 

•  1 

1 

•  1 

1 

li 

1 

CM  I  ( 

d 

1 

•  1 

1 

•  1 

1 

t» 

\ 

Ts»x  10^01%) 

5  -J 

0  ! 
X  10^  (m%) 

5  -£ 

qsgXl07(m%) 

u 

- 1 - 

1 

> 

1 

1 

1 

-0.05 

1 

1 

1 

1 

1 

1 

S-°-^ 

1 

•  ' 

1 

1 

1 

•  * 

1 

i 

1 

*• 

1 

-0.15 

•  1 

1 

•  \ 

1 

u 

1 

1 

o 

1  ro 

•  1 

1 

•  1 

l« 

1 

5  0,  „  ! 

Qb*  X  (m%) 

5  -i 

;  0,  ,  ! 
qb„x10^  (m^/s) 

5  -J 

>  0,  „  £ 
q5gXl0^(m%) 

Figure  6.8:  Variation  of  Suspended  Load  and  Bed  Load  Quantities  Over  Depth  at 
X  =  7.35  m,  Test  2  (Energetics  Model) 


Table  6.7:  Total  Sediment  Loads  With  Energetics  Model,  Test  1 
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Table  6.8:  Total  Sediment  Loads  With  Energetics  Model,  Test  2 
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CONCLUSION 


Two  irregular  wave  tests  were  conducted  on  a  movable  sand  beach  after  al¬ 
lowing  the  beach  to  shape  itself  into  equilibrium  profiles.  The  first  objective  of  these 
tests  was  to  supplement  the  results  of  Herrman  et  al.  (1997),  further  investigating 
the  cross-shore  variations  of  the  probability  distributions  and  statistics  of  the  free 
surface  elevation  and  horizontal  velocity  in  the  shoaling,  surf  and  swash  zones.  In 
addition,  velocity  statistics  were  analyzed  over  depth  at  two  cross-shore  locations. 
The  second  objective,  only  partially  achieved,  was  the  synchronization  of  velocity 
and  free  surface  time  series  for  the  separation  of  incident  and  reflected  waves.  The 
final  objective  involved  the  application  of  the  sheet  flow  and  energetics  sediment 
transport  models  to  the  equilibrium  profiles  obtained  in  this  study,  in  an  effort  to 
identify  shortcomings  of  the  models. 

The  exponential  gamma  distribution  was  again  adopted  in  this  study,  which 
examined  its  application  to  more  realistic  beach  profiles.  The  distribution  with 
the  measured  mean,  standard  deviation,  and  skewness  has  again  been  shown  to  be 
capable  of  describing  all  the  measured  probability  distributions,  though  uncertainties 
in  measurements  of  bottom  elevation  make  verification  more  difficult  in  very  shallow 
water  with  the  mean  water  depth  h  <  1  cm.  It  was  assumed  that  the  lower  limit 
of  the  exponential  distribution  of  the  free  surface  elevation,  7/,  was  imposed  by  the 
bottom  elevation  in  the  swash  zone.  This  assumption  implies  that  /frms/h  =  V^, 
where  h  =  mean  water  depth,  and  Hnas  =  root-mean-square  wave  height,  defined 
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as  i/rms  =  V^CT  with  a  =  standard  deviation  of  t].  The  measured  values  of  ifrms/^ 
in  this  experiment  increased  in  the  swash  zone  and  approached  an  upper  limit  of 
approximately  three.  With  some  irregularity  in  Test  1,  the  measured  skewness,  s,  of 
7/  generally  increased  landward,  then  decreased  as  waves  approached  the  still  water 
shoreline,  and  finally  increased  again  in  the  swash  zone  to  an  upper  limit  near  the 
predicted  s  =  2.  The  trend  in  measured  kurtosis,  K,  of  rj  was  predicted  by  the 
relation  between  K  and  s  derived  from  the  exponential  gamma  distribution,  but 
actual  data  may  require  an  empirical  fit. 

The  undertow,  u,  and  the  standard  deviation,  (T„,  of  the  horizontal  velocity, 
u,  were  again  measured  near  the  approximate  mid-depth  below  SWL  in  the  shoaling 
and  surf  zones.  In  addition,  these  velocity  statistics  were  measured  in  two  vertical 
profiles.  The  undertow  and  standard  deviation  were  expressed  in  terms  of  h  and 
Hrms  using  linear  progressive  long- wave  theory.  These  simple  relationships,  together 
with  the  measured  values  of  h  and  Hrms,  were  again  found  to  predict  the  measured 
cross-shore  variations  of  u  and  fairly  accurately  in  deeper  water.  Overprediction 
of  these  values  in  shallower  water  shoreward  of  the  bar  was  likely  to  have  resulted 
from  a  roller  effect  due  to  water  carried  onto  the  terrace  by  breaking  waves.  The 
velocity  profiles  showed  little  variation  in  u  or  cr„  over  the  depth  for  these  irregular 
wave  tests.  The  skewness,  Sy,,  of  u  was  again  found  to  be  consistently  smaller  than 
the  corresponding  skewness,  s,  of  rj,  indicating  the  vertical  decrease  of  nonlinearity. 
The  vertical  variation  of  s„  was  relatively  minimal  below  the  wave  trough  level 
and  above  the  bottom  boundary  layer.  The  kurtosis,  Ky,  of  u  exhibited  a  similar 
cross-shore  trend  to  Sy  and  also  appeared  to  remain  fairly  constant  over  depth. 

As  part  of  a  plan  to  examine  incident  and  reflected  wave  spectra  in  shallower 
water,  simultaneous  measurements  were  made  of  the  free  surface  elevation  and  the 
horizontal  velocity,  and  the  combined  incident  and  reflected  spectra  were  calculated 
at  17  cross-shore  locations.  Because  of  slightly  different  timing  rates  on  computers 


measuring  velocity  and  free-surface  time  series,  sychronization  of  the  time  series 
became  much  more  difficult  and  uncertain.  Furthermore,  the  combined  wave  spectra 
showed  little  evidence  of  long  standing  waves  generated  in  the  shallow  region  of  the 
profile,  indicating  that  the  terraced  and  barred  beach  profiles  were  not  caused  by 
low-frequency  standing  waves. 

Finally,  the  sheet  flow  and  energetics  sediment  transport  models  were  applied 
to  the  equilibrium  sand  profiles.  The  energetics  model,  which  attempts  to  account 
for  effects  of  wave  asymmetry,  undertow,  and  bottom  slope,  predicted  smaller  profile 
change  than  the  sheet  flow  model,  which  did  not  include  the  effects  of  wave  breaking 
and  undertow.  The  sheet  flow  model  is  currently  capable  of  predicting  only  onshore 
transport  outside  the  surf  zone  and  may  have  been  applicable  only  to  a  portion 
of  the  profiles  measured  in  these  tests.  Both  models  predicted  significant  changes 
in  profile  shape,  especially  around  the  region  of  wave  breaking,  and  neither  model 
was  able  to  predict  correctly  a  stable  equilibrium  "profile  using  the  velocity  and  free 
surface  data,  except  that  the  energetics  model  was  within  measurement  errors  for 
test  2. 

In  summary,  the  statistical  results  of  these  tests  indicate  that  further  devel¬ 
opment  of  a  time-averaged  model  incorporating  the  exponential  gamma  distribution 
will  be  worthwhile,  although  such  a  model  would  be  best  verified  by  detailed  field 
data.  The  model  will  need  to  predict  the  cross-shore  variations  of  h,  cr,  s,  and 
Although  they  effectively  describe  u  and  cr^  in  deeper  water,  the  simple  linear  rela¬ 
tionships  discussed  in  this  study  are  not  sufficient  for  prediction  of  these  quantities 
in  shallower  water  or  for  prediction  of  and  Ku  in  any  depth. 

The  failure  of  state-of-the-art  sediment  transport  models  to  predict  the  equi¬ 
librium  profiles  illustrates  the  strict  limitations  on  such  models  at  the  present  time. 
The  sheet  flow  model  may  be  improved  by  the  inclusion  of  additional  terms  to 


account  for  the  offshore  transport  effects  of  undertow.  The  energetics  model  over¬ 
predicted  profile  change  in  mild  wave  conditions  both  for  Thornton  et  al.  (1996) 
and  for  test  1;  it,  too,  will  need  to  describe  the  suspended  load  in  a  physically  more 
realistic  manner. 
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